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Abstract. We discuss from scratch the classical structure of Dirac spinors on an arbitrary globally 
hyperbolic, Lorentzian spacetime, their formulation as a locally covariant quantum field theory, and the 
associated notion of a Hadamard state. Eventually, we develop the notion of Wick polynomials for 
spinor fields, and we employ the latter to construct a covariantly conserved stress-energy tensor suited 
for back-reaction computations. We shall explicitly calculate its trace anomaly in particular. 
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1 Introduction 

In the last few decades we witnessed an amazing leap forward in our understanding of the 
formulation of quantum field theory on curved backgrounds thanks to the efforts of many research 
groups which often tackled this topic by means of the algebraic formalism. A careful analysis of 
all the related achievements, though tempting, would require a review on its own and, instead, 
we shall restrict ourselves to briefly mentioning the results of a recent manuscript, which has 
prompted our interest towards the topic discussed in this paper. To wit, in |DFP08j . it was 
shown that, in the framework of semiclassical Einstein's equations, it is possible to construct 
explicit solutions for a homogeneous and isotropic Friedmann-Robertson- Walker spacetime with 
flat spatial sections, where the assumption has been made that the matter content is described 
by a suitably quantised free massive scalar field. In order to prove this result, two ingredients 
have played a key role, namely, Hadamard states, as the natural candidates for a ground state, 
and the quantum behaviour of the regularised stress-energy tensor T^^. Particularly, as end 
point of the above cited paper, a late times stable de Sitter solution has been displayed and, 
hence, an effective cosmological constant has arisen without inserting it from the very beginning 
as an input datum. If one tries to seek the origin of this genuine quantum effect, one can realise 
that it is ultimately rooted in the so-called trace anomaly, i.e., in a few words, the expectation 
value of the regularised trace is not vanishing on Hadamard states, even though this is the case 
at a classical level. 

Although interesting, the derivation of the aforementioned result rises naturally the question 
about its robustness since one could wonder if this behaviour is a feature pertaining only to 
scalar fields or if it holds true for any kind of matter constituent. It thus seems advisable to 
try to apply the same scheme of reasoning in the context of a free Dirac field. To our utmost 
surprise, we have realised that the accomplishment of this goal has not been as easy as one might 
a priori believe since, just by a quick scan of the available literature, it is manifest how spinor 
fields play a somehow ancillary role in the arena of algebraic quantum field theory. In fact, 
there are only few mathematically sound results and many tools and concepts, which have been 
thoroughly discussed for scalar fields, have barely been scratched for spinors. As an example, let 
us recall that we are mostly concerned with the trace anomaly for Dirac fields and it turns out 
that this quantity has indeed been already investigated in the late seventies in |Chr78[ [ChDu79] . 
though by means of the so-called De Witt-Schwinger expansion which lacks mathematical rigour 
and, therefore, we cannot call our understanding of this anomaly complete yet. 
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With this in mind, before we can tackle any specific topic in a cosmological scenario, our first 
concern must lie in the amendment of the the above mentioned problem. This will indeed be the 
main point of the paper and we shall discuss it within the framework of the algebraic approach to 
quantum field theory. To this end, we must also take into account that the scientific community 
interested in this problem might not be acquainted with the formulation of a spinor field theory 
on a curved background, which, already at a classical level, turns out to be rather different and 
more complicated than the scalar counterpart. Hence, as a starting point, we shall review the 
construction of a Dirac field in a classical framework emphasising the role of the underlying 
geometric structures which are needed in order to fully describe both the kinematically and the 
dynamically allowed field configurations. To this effect, our analysis will benefit from earlier 
works which have already dwelled on this topic and, most notably, we shall refer to the seminal 
paper of Lichnerowicz |Lic64] as well as to |Dim82l ISaOSj . 

Subsequently, we shall discuss the quantisation of a Dirac field on a curved background and, 
in this respect, one should mention that there are several possibilities at our disposal. On the one 
hand one could follow the point of view already suggested in [Dim82| . while, on the other hand, 
one could also analyse the problem from the perspective of Araki |Ar70j , whose scheme has the 
peculiarity of unifying spinors and cospinors in a single body before quantising them. This leads 
to a natural definition both of a CAR *-algebra of fields and of a subalgebra of observables, 
once we require at least that elements whose supports are spacelike separated must commute. 
Furthermore, this scheme, also at the heart of [Ver96[ lKo95l IKraOOi IHol99[ IFeVe02l ISaOS] . has 
the remarkable advantage of being well suited to recast the quantum theory of Dirac fields in 
the language of Locally Covariant Quantum Field Theory |BFV03| . as we will point out. 

In order to fully control the machinery of a quantum field, the scalar field scenario already 
thought us that one has to understand well which algebraic states one should use and, to this 
avail, the ones of Hadamard type are the natural choice in the context of Dirac fields, too; 
these play the role of ground states in a curved background and their ultraviolet behaviour 
closely mimics that of the Minkowski spacetime vacuum state. Consequently, the fluctuations 
of the components of the stress-energy tensor on these states are bounded, a property which 
is vital in the context of semiclassical Einstein's equations. Hadamard states for Dirac fields 
have already been discussed in [K5951 IVctMI IHo199[ HMnl ISaVenil IDaHoOGj and we will review 
their properties in detail before employing them to achieve the first of our main results, namely, 
the construction of the extended algebra of Wick polynomials in a spinor framework. To this 
end, we will follow the path paved in the scalar scenario in |BFK951 IBrFrOOl IHoWaOl] IHoWa02] 
where it has been displayed that such polynomials lie at the basis of a sound S-matrix formalism 
for interacting field theories on a globally hyperbolic curved background. 

Notwithstanding, Wick polynomials are already valuable in free field theories and we will, as 
already anticipated, indeed use them to achieve our second main result, i.e., the definition of a 
well-behaved quantum stress-energy tensor operator T^y for Dirac fields. To achieve this goal, we 
shall follow a procedure similar to that discussed in |Mo03| (see also the related work |HoWa05] ) 
for a free scalar field, i.e., we shall introduce an improved point-splitting procedure to define 
T^y evaluated on a Hadamard state. This leads to a new overall stress-energy tensor which does 
not alter classical dynamics and is ultimately conserved at the quantum level. Furthermore, as 
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a by-product of our analysis, we shall also be able to explicitly compute the expectation value of 
its trace which will agree, up to terms proportional to OR, with previously found results while 
being derived in a rigorous framework. 

2 Dirac fields: a classical overview 

Since, as we have outlined in the introduction, the aim of this paper is to provide an as much 
as possible self-contained approach to some topics related to the quantum description of Dirac 
fields in curved backgrounds, we will start with a description of Dirac spinors in a classical 
framework. Although such topic has been already discussed both from a geometrical and from 
an analytical point of view by many authors, we reckon the we should try to recall the main 
features of the classical approach in order both to facilitate the understanding of the quantum 
aspects and to fix some subtleties which ubiquitously arise in these scenarios. 

2.1 On the spin structure and related geometric entities 

Bearing in mind this overall philosophy, we shall mainly devote this subsection to the introduc- 
tion of spin structures and of the Dirac bundles, in order to characterise (co)spinors as suitable 
sections. We shall not dwell too much on the geometrical contents and for the potential readers 
who might find our approach too shallow we present our apologises and point them to [LaMi89j 
for a careful discussion of most of the forthcoming concepts and applications. 

As a starting point, let us fix that, in this paper, a spacetime is meant to be a four- 
dimensional, Hausdorff, smooth manifold endowed with a Lorentzian metric, whose signature is 
chosen as (— , -|-, -|-, +). Furthermore, since it is common wisdom to associate Dirac fields to the 
notion of spin, we need a few definitions as a first step: 

Definition 2.1. We call spin group Spin{p,q) with p, g € N the double cover of SO{p,q), 
i.e., it exists the following short exact sequence: 

{e} — >Z2 — > Spin{p, q) — > SO{p, q) — > {e} , 

where {e} stands for the trivial group, whereas Z2 = {±1} is the cyclic group of order 2. There- 
fore, any element of Spin(p,q) induces an element of SO(p,q). Such a surjective covering will be 
indicated as 11 : Spin{p,q) — > SO{p,q). 

Remark. As a consequence of the above definition, it exists an isomorphism between Spin{p, q) 
and Spin(q,p) for any possible value of both p and q. Furthermore, it is also possible to talk 
about the dimension of such classical Lie group which, per direct inspection, is 

dim {Spin{p, q)) = ^{p + q){p + q - I) ■ 

Furthermore, for all p,q > 0, the spin group has two connected components, where we denote the 
component connected to the identity as SpinQ{q,p). The latter insight entails that the scenario 
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withp = 3 and g = 1 is of great interest since, in this case, Spino{l, 3) is isomorphic to SL{2, C). 

The above definition represented only the first step toward the definition of a Dirac field 

since, in modern classical field theory, the geometric interpretation of a kinematically allowed 
configuration is that of a section of a suitable associated bundle. Within this respect, one should 
notice that, (un)fortunately, in the often analysed case of a scalar field, such nice perspective 
does not really enter the fray, whereas, in this case, such a luxury is not at our disposal, being a 
spinor intrinsically a vector- valued field. Hence, we shall now show how the notion of spin group 
can be intertwined with that of a classical field. As a first step we need further definitions. 

Definition 2.2. Given a vector bundle E = E[V, tt, M] over a Lorentzian manifold M with 
typical fibre V and projection map tt : E ^ M, we call frame F{x) over the point x the 
assignment of an ordered basis to the fibre 'k~^{x) = V , i.e., a map p : 7r~^(x), being k 

the dimension ofV and K € {M, C}. Furthermore, we require that the right action of the group 
GL{k,K.) on the fibre is both free and transitive. 

Remark. Such a definition is at the heart of the widely exploited concept of tetrads in general 
relativity and, on a practical ground, it is remarkable since it guarantees us that whenever a 
vector bundle is associated to the underlying spacetime one is free to choose a basis of such a 
space and that the rules of changing the basis are left untouched; this is encoded in the action of 
GL[k, K), i.e., each element of this group transforms a basis into another one. Furthermore, the 
request of transitivity guarantees us that it is always possible to transform any chosen basis into 
any other one, whereas demanding a free action entails that only the identity element leaves a 
chosen basis unchanged. The most notable and ubiquitously used application of such a definition 
is the tangent bundle where V = M.^, being k = dimM. In the following, we shall have this case 
in mind and hence we shall identify E as, E = TM = TM\^^, vr, M]. 

Therefore, if we call F^M ~ GL{k,'R) the set of all possible frames over a point x G M, we 
can gather all this information into a unique object: 

Definition 2.3. A frame bundle FM associated to TM is FM = Ux^xM. This is a fibre 
bundle on its own, namely FM = FM[GL{k,R),TT' , M], with GL{k,W) as typical fibre and the 
projection map it' : FM M . Whenever the base space M is an orientable and time orientable 
spacetime, hence also endowed with a (local) non- degenerate metric of signature {p,q), we can 
choose a space and time orientation, and hence reduce GL{k,R) to SOo(j),q) withp + q = k. 
We shall always consider such case in the following. 

Remark. A fibre bundle as above is an example of a principle fibre bundle, having the same 
group as both the typical fibre and the structure group. By requiring the right action of the 
group to leave the base point invariant, e.g., t^'oRck^^^^^ = tt' for FM, one can straightforwardly 
extend it to the full principle fibre bundle. This extension is independent of local trivialisations, 
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since these are related by the left action of the group and the left and right action of a group 
on itself commute. One can even switch perspectives and define the fibres and base points of a 
principle fibre bundle via the orbits of the right action. 

We are now in the position to eventually introduce the main geometric structure of this 
paper which lies at the heart of the construction and of the analysis of a Dirac (co)spinor: 

Definition 2.4. Given a manifold M with a non degenerate metric of signature {p, q), a spin 
structure is the pair {SM,p) where SM = SM[SpinQ{p, q),7r, M] is a principle fibre bundle 
over M with the identity component of the spin group as typical fibre. Moreover, p is a bundle 
homomorphism from SM to FM subject to the following conditions: 

1. p is base point preserving, such that 

Tt' O p = TT, 

2. p must be equivariant, i.e., calling i?^ and the natural right actions of SpinQ{p,q) on 
SM and of SOo{p,q) on FM respectively, we require that 

p o i?^ = o p, VA G SOo{p, q), 

where A = n(A), being H as in definition \ 2. 1[ 

Remark. A natural, apparently nai've, but ultimately rather complicated question a potential 
reader might ask, is why one should cope with such complicated geometric structures to deal with 
a somehow natural concept such as that of spin. To answer this question, let us keep in mind 
the notion of a Dirac field in the fiat scenario and try to generalise it to a curved background 
with Occam's razor as a principle. That said, in four-dimensional Minkowski spacetime, one 
can employ a standard construction, dating from the late thirties and due to Wigner, according 
to which a spin ^ field is nothing but a suitable function satisfying the Dirac equation and 
transforming under a unitary and parity invariant representation of the universal cover of the 
Poincare group induced from the j = \ representation of SU{2) (see chapter 21 of |BaRa86j 
for more details). It is remarkable to notice that, once restricted to SL{2,C), i.e., neglecting 
the translations, this representation coincides with the D^^'^^ © D^^'^^ representation of the 
aforementioned group. Although there is no translational invariance in a curved background, 
the previous definitions grants us that, by means of the spin structure, it is nonetheless possible to 
define, barring certain technical requirements, a natural notion of an 5^(2, C) -group associated 
to a spacetime. Furthermore, the very absence of translational invariance implies that such a 
definition of spin in curved backgrounds seems really the best one can do. 

Hence, though not sufficient to determine full dynamics of a Dirac spinor, our philosophy 
will be to seek to characterise the kinematically allowed configurations of such a field by means 
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of the mentioned £)(2'°) D^^'^) representation, while remembering at the same time that, 
classically, fields should be understood as sections of a suitable vector bundle. To combine these 
two concepts, we proceed in the following way: 

Definition 2.5. We call Dirac bundle of a four dimensional Lorentzian spacetime M with 
respect to the representation T = d(2'°) © D^^'^) of SL(2,C) on the associated bundle 
DM = SM Xt C^. This is the set of equivalence classes [{p,z)], where p G SM, z & and 
equivalence is defined out of the relation 

{pi,zi) ~ (^2,2:2), 

if and only if it exists an element A of SL{2,C) such that = p2 and T{A~^)zi = zi, 

where La denotes the left action of A on SM . The global structure of DM is that of a fibre 
bundle over M with typical fibre C^, and the projection map ttd is traded from the one of SM, 
namely, V [{p, z)] G DM, it holds 

t^d[{p,z')] =tt{p). 

Furthermore, if we endow with the standard nan degenerate internal product, we can construct 
the dual Dirac bundle D*M as the C^* -bundle associated to SM requiring that and 

{La{pi), z*T{A)) are equivalent, where * denotes the adjoint with respect to the inner product 
on and elements of C^* are understood as row vectors. Consequently, the dual pairing of 
and C^* extends in a well-defined way to a fibrewise dual pairing of DM and D*M . 

Definition 2.6. Let E be an arbitrary vector bundle over M. With S-{E) = C°°{M,E), we 
denote the space of smooth sections of E, endowed with the topology induced by the family 
of seminorms 

\\f\\k,c = sup{\f'^'^\x)\,xeC}, fee.{E), 

where C is an arbitrary com,pact set and ■ ^''^ denotes a derivative of k-th order. Furthermore, 
we introduce the space of smooth sections with compact support T)(E) = Cq-'{M,E), 
equipped with the topology induced by the family of seminorms 

\\f\\k = sup{\f^^\x)\}, feD{E). 

In the case E = DM, E = D*M, we can define a global pairing of8.{D*M) and D{DM) 
or D{D*M) and 8,{DM) by integrating the local pairing induced by the inner product on C^, 
e.g., 

if, 9) = j dfi{x)f{x) {9{x)) 

M 

for all f e E{D*M), g eD{DM). 

We are finally in the position to define the key object of our analysis: 
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Definition 2.7. A Dirac spinor is a smooth global section of the Dirac bundle, i.e., ip € 
£.(DM) or, equivalently, if we consider a local open neighbourhood U of any point x G M , ^ is 
(dijjeomorphic to) a A-vector field tpUj i-^-y '4'\u ~ V'c/ '■ U ^ C^, since the bundle trivialises as 
DM\ij ~ {/ X C^. Analogously, we call Dirac cospinor a smooth global section of the dual Dirac 
bundle, namely, 'ip' G E{D*M) or, in a local neighbourhood U, ip'lu ~ tpij '■ U C^* ~ C^. 

We would like to stress that our definition of Dirac spinor fields does not include any equation 
of motion. 

It is quite safe to admit that the introduced geometric objects are rather technically compli- 
cated, and one might wonder if the class of spacetimes admitting them is not rather restricted. 
Particularly, from a physical point of view, one is interested in introducing spinors and fields 
in general as global objects. There is thus a compelling need to understand which is really the 
set of backgrounds we can work with, and the answer is surprising and reassuring at the same 
time. In fact, the following theorem can be proved (see [BoHiSSj for the original proof and also 
|Ge68j for a more physically oriented analysis and proof): 

Theorem 2.1. A manifold M admits a spin structure if and only if it has a vanishing second 
Stief el- Whitney class or, in other words, if the second de Rhan^ cohomology class i?^(M, Z2) 
is trivial. Furthermore, if the manifold is four- dimensional, space and time-oriented, as well as 
globally hyperbolic, this requirement is automatically satisfied. 

This theorem is rather useful because globally hyperbolic spacetimes are the most interesting 
and natural class of manifolds whenever one deals with both classical and quantum field theory 
on curved backgrounds. As a matter of fact, each of these spacetimes can be foliated as S x M, 
being S a smooth Cauchy surface [BeSaOSj . Therefore, on these backgrounds one can state 
precisely the notion of initial value problem for the equations of motion, hence determining 
the classical dynamically allowed configurations of a field as the solution of a certain partial 
differential equation. 

The property of globally hyperbolic, four-dimensional spacetimes M which guarantees the 
existence of a spin structure is their parallelisability, i.e., the fact that there always exists a 
global orthogonal frame on them. Consequently, FM is a trivial bundle in that case, and this 
property extends to SM, TM, T*M, DM, and D*M as well. We are thus in the position to 
introduce global frames of the latter four bundles. 

1. Employing a global section E of SM, we can define a spin frame {Ea}a=i---4, ^-c, a 
set of four global sections of DM as Ea{x) = [(E{x), za)], being za the standard basis 
of C^. Hence, any Dirac spinor ip can be decomposed as 'ip{x) = ^l^^{x)EA{x) where now 

2. A dual spin frame {E^}b=i,— ,4, i.e., a set of four global sections of D*A'I can then be 

^In the most general framework one should consider the second Cech cohomology class, but this coincides with 
the de Rham one for differentiable manifolds. In this paper we deal solely with spacetimes which, according to 
the definition stated at the beginning of the section, are differentiable. 
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automatically constructed out of the frame requiring Ej^{E^) = 5^. From now on capital 
letter indices will refer to quantities expressed with respect to these bases. 

3. Exploiting definition \'1A\ we can project E to FM, hence obtaining a global section e = 
p o E oi EM. Employing this, we can define a Lorentz frame {ea}a=o,.-- ,3, i-e., a set 
of four global sections of TM, by realising that TM can be understood as the M^-bundle 
associated to EM. Such a fibrewise basis is orthonormal in the sense that g{ea, et) = rjab, 
where r]ab is the Minkowski metric, and it is often referred to as the non-holonomic basis 
of the base manifold; this is in sharp contrast with the standard (holonomic) coordinate 
basis which is related to the basis by means of a basis change or, in other words, by 
the matrix Ca denoting the coefficients of in their expansion with respect to d^. 

4. Analogous to the spin case, one can now straightforwardly define a dual Lorentz frame 
{e^}f,=o,... ,3 constructed out of Ca as Ca(6^) — '^a' From now on, lower-case Roman let- 
ters will always refer to quantities expressed with respect to the non-holonomic basis, 
whereas lower-case Greek ones will indicate those with respect to the holonomic basis. 
Non-holonomic indices will be "raised" and "lowered" with rjab, while the same operations 
will be performed on the holonomic ones using g^y = g[d^,di,). 

The most notable consequence of this new data is that we can decompose evert spinor-tensor 

f G £(TM ■ ■ ■ (g) TM ® T*M (g) ■ ■ ■ (g) T*M (g) DM • • • ® DM ® D*M • • • » D*M) 

i i k I 

as follows 

/ = fhl'l-k'-Bi'''^^ ®---®ea,®e^' (^---^e^^ ®Ea^®---®Ea^®E^^ ®---®E^\ 

One could in principle certainly choose a different global section E' of SM and thus obtain 
different spin and Lorentz frames which are related to the previous ones by local spin and 
Lorentz transformations. On the level of coefficients, such a change of frames results in 

where A G r(SL(2,C)), whereas A = n(A) G 50o(3,l). 
2.2 On the dynamics of a classical Dirac field 

Since we have by now assured ourselves of the existence and well-posedness of the global structure 
of Dirac fields in a globally hyperbolic time and space-oriented manifold, we shall next proceed 
to introduce the natural evolution operator out of which one can describe the classical dynamical 
content of our theory. It is imperative to stress a sharp difference between the previous and 
this subsection; in the preceding discussion, all the introduced geometric structures have been 
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somehow natural and intrinsic, i.e., no special choice has been performed with the due exception 
of theL>(5'0)eD(0'3) representation to define the Dirac bundle. Conversely, in the forthcoming 
discussion, some arbitrariness appears and we shall try to emphasise it to a potential reader, 
since we have to keep track of it to ensure that it does not play a distinguished role in the 
forthcoming discussion of the quantum fields. To wit, we are referring to the definition of the 
so-called 7-matrices. To obtain them, we can proceed as follows: 

Definition 2.8. Given M^'"^ endowed with the metric rj of signature {p, q), we call Clifford al- 
gebra Cl{p, q) ofMP''^ the real associative algebra generated by the identity I and an orthonormal 
basis ofM^'^ whose elements 7^ with a = 1, ...,p + q are subject to the relations 

{7a, 7b} = 7a7fe + 7fe7a = '^Vabl- (1) 

Particularly, if p = 3 and q = 1 or vice-versa, one can refer to (7/(3,1) or C/(l,3) as Dirac- 
Clifford algebra. 

Remark. It is a direct consequence of this definition that a basis for the Clifford algebra is given 
by the identity and by all products 7ai----7a„ with oi < ... < a„ and n < p + q, which entails 
that dim{Cl{p, q)) = 2^+'^. As a further important datum, we wish to underline that Cl{p,q) is 
a Z2-graded algebra; this arises if we introduce the automorphism a : Cl{p,q) — > Cl{p,q) such 
that Q!(7a) = —7a for all possible a. Since coincides with the identity map, we can always 
decompose: 

Cl{p,q)=Cl\p,q)(BCl\p,q), 

where Cr{p,q) = {a G Cl{p,q) \ a{a) = (— )*a}. By direct inspection, one can realise that 
Cl^{p, q) is the subalgebra of the full Clifford algebra generated by products of even numbers of 

7a- 

The Dirac-Clifford algebra enjoys many relevant properties of great interest for our dis- 
cussion. As a first step, using the periodicity theorem 4.1 in |LaMi89j . one can prove per 
direct inspection that CZ(l,3)c = C/(l,3) (8) C is isomorphic to the algebra M(4, C) of 
4x4 complex matrices. This entails that it is natural to seek for a complex representation 
T : CZ(l,3)c Hom{C^ and, according to theorem 5.7, still in [LaMi89] . there is only 
one of these which is, up to equivalence, irreducible. Its matrix form can be described by the 
matrices 7^, which we choose in such a way that (7o)'^b ~ ~'^0B whereas (7^)'^^ = 7^ for 
a = 1, 2, 3, i.e., 

where a = 1,...,3 and a a is a Pauli matrix, whereas /„ denotes the n x n- identity matrix. 
Furthermore, these matrices, independently of the choice in ([2]), are the so-called ^-matrices and 
they can always be interpreted as the coefficients of a global tensor 7 E 8.{T*M ® DM ^ D*M) 
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admitting the following expansion: 



7 = lase" ®Ea0E 



(3) 



This identity yields that, if we use the tetrad coefficients to define 7^ = 7ae^, we recover the 
standard anticommutation relations for the 7-matrices in curved backgrounds 



This is the net effect of changing from a non-holonomic to a holonomic basis. 

Let us stress at this point that we have chosen a specific representation of the Dirac-Clifford 
algebra out of the possible equivalent ones and indeed this is the arising arbitrariness announced 
at the beginning of this subsection. The choice of different representations of the Dirac-Clifford 
algebra can be shown to lead to quantum field theories which are equivalent up to gauge trans- 
formations [Sa08j and, moreover, if one restricts to the quantum observables, the choice of a 
representation becomes even irrelevant, as we will discuss in the next section. 

A second interesting property of the Clifford algebra arises from the realisation that it con- 
tains the spin group (see §2 of [LaMi89j ) and, most notably, Spin{3,l) C C/o(3, 1). It is thus 
natural to wonder how the above introduced representation yielding the 7-matrices restricts to 
the spin group, and the answer to this query comes from proposition 5.15, still in [LaMi89] . 
which guarantees us that each irreducible complex representation of the Clifford algebra on a 
vector space can be restricted to the sum of two inequivalent irreducible representations of the 
spin group. In the case under analysis, one can by direct inspection realise that the restriction 
of T to Spin{3, 1) yields the sum of two non-equivalent irreducible representations, which, on 
51/(2, C) ~ Spino{3, 1), coincide with the previously mentioned D^^^'^) Q) D^^'^^ representation. 
According to our analysis, for any vector field v G 8,(TM), we can meaningfully introduce 

V = V^la, 

which is an element of E,{DM 0D*M), such that its coefficients v^^ form a 4 x 4 complex matrix. 
Notice that, from now on, we will not specify explicitly when we shall deal with abstract Dirac- 
Clifford algebra elements 7a or with their matrix representations. It is understood that, whenever 
we either contract 7-matrices with a vector field or these matrices are applied to a vector in C^, 
we refer to the latter case. 

The last ingredient we need to specify the dynamics of Dirac fields is a parallel transport 
on the Dirac bundle. The grand strategy is rather simple, namely, we introduce the standard 
metric connection, interpret it on the frame bundle and eventually lift it to both the spin and 
the Dirac bundle: 

Definition 2.9. Let lo : FM — > T*FM o(3, 1) denote the connection 1-form of the unique 
Levi-Civita connection on FM . It induces the standard Levi-Civita connection on TM (and 
vice versa) which can be expressed as the covariant derivative 



(4) 



V : E{TM) &{TM (g) T*M), Ve^ = K^eb 
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where the connection coefficients T^^ of the Levi-Civita connection are specified by 

rL = ((^ ° e) [e* o Ca] Cc) , 

with : T*M T*FM denoting the push-forward of e in the sense of cotangent vectors. 

The pull-back Q = (dll)"-^ o p* o uj of uj to SM, with dU : 5/(2, C) — > o(3, 1) denoting the 
derivative of the covering H at the identity, defines the spin connection, which by the definition 
of DM as a bundle associated to SM can be specified as a covariant derivative 

V : E{DM) ^ E{DM ®T*M), VEa = (t^a^^ Eb, 

where the spin connection coefficients are given by 

a^A = E''{{noE)[E,oea]EA) (5) 

and E^ : T*M T*SM denotes the push-forward of E in the sense of cotangent vectors. 

We will soon prove that the coefficients a^js^ can be expressed in a simple way by means of both 
the coefficients and the 7-matrices. Before we tackle such task, let us stress that the covariant 
derivatives V from definition 12.91 can be straightforwardly extended to any spinor-tensor 

/ G £(TM • • • TM ® T*M ® ■ ■ ■ T*M (g) DM • • • ® DM ® D*M • • • D*M) 
by defining 



V : £(rM 



>TM(E)T*M' 



>T*M(S)DM 



DM (g) D*M ( 



D*M) 



E{T*M(^TM> 



TM(S>T*Mi 



k I 

T*M ® DM • • • DM D*M • • • (g) D*M). 

' V ' V ' 

k I 

(6) 



At a level of components and, for notational simplicity, only in the case i = j = I = k 
reads 

V/ = e'Vc (fl^iea ®e^®EA®E^ 



1, 



Q faA T^d faA , -pa fdA „A faC i _C f 
UcJbB ~ ^ cbJdB "T cdJbB ~ ^cCJbB ^cB// 



bC 



' e'' Ea E^ . 



As promised, we shall now give an explicit expression for the (j-coefficients. Although a 
demonstration of the next lemma is already present in |Lic64j (see also chapter 13 of |Wa84j ) . 
we shall prove it again due to ubiquitous sign subtleties arising from the choice of the metric 
signature and possibly leading to confusions when comparing with the literature, |Dim82j in 
particular. 
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Lemma 2.1. The connection coefficients of the spin connection can be expressed as 

l-pft B dC (^\ 
^aA = -^adihCl A- V) 



Proof. The strategy of the proof is the following: we first derive an explicit expression for the 
double covering homomorphism 11 : Spino(3, 1) ~ 5*^(2, C) S0o{3, 1). From this we obtain 
an expression for its derivative at the identity dll : 5/(2, C) o(3,l), which, inserted in ([5]), 
yields the wished- for result. 

Let us thus recall that Spino{3, 1) can be understood as a subgroup of CZ(3, 1). This entails 
that, for any A G SL{2, C), we can define the adjoint action 

Being SL{2, C) a finite cover of 5*00(3, 1), applying an induction-reduction mechanism, we know 
it exists a representation T of 50o(3, 1) such that, setting A = n(A), we have 



T{K)^a = Ad~A^a)- 



(8) 



Since we are dealing with Lie groups, we recall that all finite-dimensional representations are 
matrix representations, i.e., we can simply write T(A)7a = A'^^^jh- Hence, A7aA~"'^ = A^^-fb, 
where we notice the invariance of the left hand side under the Z2-action sending A — > itA, as 
one could have expected, being SL{2,C) the double cover of 5*00(3, 1). 

Let us now take an arbitrary differentiable path t i-^ A{t) in 5L(2,C) whose projection on 
50o(3, 1) is differentiable; the following identity holds 

If we derive with respect to t, an identity between algebra representations arises, namely, 

dA{t) ~, , 1 ~, , dA(t)-i /dA{t)\'' 



dA{t) 
dt 



7aA(t)-i-A(t)7«A(t)- 



dt 

_^dA{t] 
dt 



Ait)-' 



lb, 



(9) 



where we have exploited that the derivation of A{t)A(t) ' = I yields 



dA{t) 



A{t)-' = -A{t) 



dA{t)-' 



dt ' ' dt 

Applying the adjoint action of A{t)~' to we end up with 



(dmy 

\ dt J ^ 



lb 
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We use ([8]) as well as the basic property of a representation, namely, Ad{A ^) = Ad{A) ^, to 
derive 

If we call A = A{t)^^^^^^^ and /i = A(t)~^^^^, it holds — 7aA = ^jf' alb- Let us multiply this 
identity on the right with 7^ and, if we bear in mind that 7"7a = rj°'^^alb = v'^^Vab = 4, we end 
up with 

4A - 7aA7'^ = llabl^l'. (10) 

Taking into account the antisymmetry of ^ € o(3, 1) and the identity 7"7'''7'^^a = 0, a possible 
solution of (fTOj) is 

A = ^/^afe7 7 • 

The uniqueness of this solution is not guaranteed, since, being the left hand side of (jlOp linear 
in A, we have merely found a particular solution, and we are free to add any further solution 
of the homogeneous counterpart, i.e., any A such that 4A — 7aA7" = 0. This implies that the 
commutator [7a, A] = 0, and Schur's lemma (see th. 4.26 in |Ha03] ) entails that A = kl, being 
/ the identity. 

The value of k can be unambiguously determined if we notice that, according to our previous 
construction, A is an element in the algebra of SL{2, C), which consists of matrices with vanishing 
trace. Therefore, if we impose Tr(A) = Tr{kl) = 0, the only possibility is A; = 0. 

Since the differentiable path chosen in the proof has been arbitrary, we have now proven the 
explicit form of dU in terms of its inverse, namely, 

(dU)-^ : 0(3,1) ^s/(2,C), {dUr\fiab) = ll^ab^l' y/^ab S o(3,l). 

Remembering the definition of 17, recalling e = po E, and inserting the expression of (dH)^^ 
into ([5]), we finally obtain 

1 rb ^dc 
(^aA = ^iad7fec7 A- 

□ 

Lemma 2.2. The Dirac- Clifford ^-matrices are covariantly constant, i.e., V7 = 0. 

Proof. This is a straightforward calculation, once a subtlety has been clarified: the matrices 
7 are constructed as the matrix form of an irreducible representation of the Clifford algebra 
and then subsequently glued to each point of the underlying base manifold M via ([3]). This 
prescription entails that da^b = for any a,b = 0, 3. 
We can now compute 

V7 = eVfe (7^56'^ (^Ea® , 
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which, exploiting definition 12.91 becomes 

V7 = [abT^se" ^Ea(S)E''+ 7^^ {^1^' ® 
-a^Ae" (^Ec®E^ + aj^cc" (^Ea®E^')]. 

If we apply formula ([6]) to sections of T*M ®DM ®D*M and consider that the matrix elements 
are constant functions, the lemma is proved out of direct substitution. □ 

We are ready, at last, to describe the dynamically allowed configurations of a spinor field in a 
curved background. Particularly, if we now understand multiplication with the 7-matrices to 
act from the left on spinors and from the right on cospinors, we shall call dynamically allowed 
a Dirac (co)spinor ip^'^ which satisfies the Dirac equatioril 

= (_y + m) V = 0, (11) 
D'^' = (y + m) V'' = 0, (12) 

where e 8,{DM), whereas i>' G Z{D*M). 

2.3 The Cauchy problem and the fundamental solutions 

In this subsection we shall discuss the classical initial value problem for Dirac fields. As it is 
well-known from fiat spacetimes, a solution of the Dirac equation is usually related to a solution 
of a hyperbolic differential equation. 

Lemma 2.3. Every solution ijj of ()lll) is also a solution of the spinorial Klein- Gordon 
equation 

= (vaV" - I - V = 0, (13) 

where R is the scalar curvature of {M,g). A similar statement holds for each cospinor ip' , 
solution of (fT2]) . 

Proof. Let us take a solution of (jlip and let us multiply it with D' . We end up with 

D'Dip = (y + m)(-y + m)tp = {-fY + m'^)^p = 0, 
which means that we need to prove that y^ = VaV^ — -j. To this avail, let us write 

= 7"V„7'Vb = ri'^a^b = 7^"7''^VaV6 + 7t"7'^ VaVfe, 

^ In Minkowski spacetime, the Dirac equation is most notable for having a purely imaginary coefficient i in 
front of the Dirac operator. Here, such a number does not appear and the underlying reason is rooted in the 
employed convention for both the metric signature and the sign of the defining anticommutation relations for the 
Clifford algebra. To wit, i is absorbed in the definition of the 7-matrices, as stated in (|2}. 



15 



where in the second equahty we have exploited lemma [2^21 whereas in the third one we split the 
expression in its symmetric and antisymmetric part. Since 7(^7'') = ^{7^,7''} is equal to the 
metric times the identity, it holds 

y2 ^ y^ya ^ ^[a^^'l v,Vb = V.V^ + 7"7'V[,V,] = VaV^ + ^7"7'Ca6, (14) 

with Cab denoting the curvature tensor of the spin connection. 

Let us briefly state some properties of Cab, which are examined in more detail in appendix 
lA.ll firstly, from lemma [2?T] one can infer that 

Cah = ■^Rabcdl'^Y ■ 

Employing the Clifford relations and the symmetry properties of the Riemann tensor Rabcd, it 
is straightforward to show that 

1 /? 

7 7 (-^afe = -7 7 ^ab = --jl 7 i^ab = 

Inserting this into p4|) . we finally obtain 
and thus 

D'D = DD' = -P= - VaV" + ^ + m^. 

□ 

Remarks. One should notice that the above lemma can be seen as a particular application of 
Weizenbock's formula [LaMi89] . 

It is furthermore worth noting that VqV" in the above expression is not diagonal in the 
spinor indices, and thus not □ times the identity. Its principal symbol g^^kfj^k^, however, is 
indeed diagonal and even of metric type. 

Let us also stress that, in sharp contrast to the scalar Klein-Gordon equation in four space- 
time dimensions, there is no freedom to select a coupling for Dirac fields to the scalar curvature 
since this is universally fixed to \. Furthermore, for Dirac fields, such factor corresponds to the 
conformal coupling whereas, for the scalar case, the latter is g. 

The introduction of the Dirac operator and the discussion of its main properties allow us to 
state and to prove the main theorem related to the classical dynamical behaviour of spinors on 
curved backgrounds: 

Theorem 2.2. Let {M,g) be a four-dimensional globally hyperbolic oriented and time oriented 
spacetime, such that M is diffeomorphic to S x M, with S a three-dimensional Riemannian 
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manifold. Let us furthermore call l : ^ M a smooth embedding of S into M , such that 

is a Cauchy surface of M . If we refer to DT, as the bundle constructed out of the pull-back of 

DM via i, then the following Cauchy problem admits a unique solution 



".1 = " . (15) 



where ip £ e{DM) and ipo G D(L>S). 



Proof. For notational simplicity we will omit the embedding l and just write ^js instead of 
for any ip G ^{DM). 

According to lemma 12.31 each solution of the Dirac equation also solves a spinorial coun- 
terpart of the Klein-Gordon equation, namely, Pip = (□ — j — m^) ^ = where P has metric 
principal symbol and is thus a normally hyperbolic differential operator. 

Hence, we can invoke the results on hyperbolic partial differential equations (see |BGP07j 
and theorem 3.2.11 in particular), which guarantee us that any Cauchy problem for compactly 
supported initial data for the the considered partial differential equation admits a unique smooth 
solution supported in the causal past and future of the initial datum. The only problem left is to 
give a prescription on how to switch from a Cauchy problem for the Dirac equation, hence with 
only one given initial datum, to one for a second order hyperbolic partial differential equation, 
where two data on the Cauchy surface have to be prescribed. To solve this dilemma, let us not 
deal with (llSp . but with the following system: 

f (□ - f - m2) u = 

I n|s = 0, -|h|s = ^Vo ' 

where u G E-{DM) and n denotes the vector field normal to S such that rj"''^nani, = — 1. As stated 
before, such a system admits a unique global smooth solution H and, hence, let us introduce 
the smooth section tp = D'u. It is immediate to see that ip is a solution of the Dirac equation 
= 0. Since ip is smooth, it can be directly evaluated on the Cauchy surface T, where 
ipls = D'u\-£ = {y + m) u\y;. Inserting the initial condition for u, the second term vanishes, 
whereas the first reads 

= -l^na-^^ni^ipQ = -r]°-''nanbipo = ipQ. 

E 

Hence, the section ip constitutes a unique solution of the Cauchy problem in the thesis of the 
theorem. □ 

Since our aim is to deal with quantum field theory over curved backgrounds in the long run, 
it is more useful to prove not only the existence and uniqueness of the solution of a Cauchy 
problem, but also the existence of the so-called fundamental solution, which is the last relevant 
theorem of the classical theory that we shall need: 



dn 
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Theorem 2.3. The Dirac operator admits unique advanced ) and retarded ) fundamental 
solutions, i.e., continuous linear maps : D{DM) — > 8.{DM) fulfilling DS^ = I = S^D. 
These maps are determined by their support properties 

supp{S^f) C J^{supp /), V/ G D{DM), 

with J^{U) denoting the causal future /past of the set U. Similarly, there exist unique advanced 
and retarded fundamental solutions : T){D*M) £{D*M) of D' . 

Proof. The strategy wih be similar to the one employed in the proof of the existence of the 
solution of the Cauchy problem. Thus, let us start from P = —D'D = —DD'; since this is a 
normally hyperbolic differential operator, we already know (see |BGP07t IHo94j ) that it exists a 
unique advanced - say E'^ - and a unique retarded - say E~ - fundamental solution for P on 
ViDM). Hence, for any / G Ti{DM), we know that 

PE^ = 1 = E^P and supp{E^f) C J^{supp /). 

We can now define = —D'E^ which, per direct inspection, satisfies DS^ = I and is further- 
more continuous and has the correct support properties, since the application of D' preserves 
these features. In the same way, we can construct advanced and retarded right fundamental 
solutions for the dual Dirac operator as 5^ = —DE^, with Ef^ being the fundamental solutions 
of P on V{D*M). 

The next step consists of proving that the right fundamental solution is "also a left one" and 
we only show this for S*^, since the proof for is analogous. Consider any h € CD(Z)*M) and 
any / G V{DM). Since Df G D(L'M), we end up with 

{h,S^Df) = {D'SfKS^Df) = {S^h,DS^Df) = 
= {S^h,Df) = {D'S^h,f) = {hJ), 

where all expressions are well-defined since suppS'^/fl supp5J/i is compact due to global hy- 
perbolicity of M. 

It remains to be shown that the fundamental solutions are unique and again we only prove 
this for 5=^ here. To this end, let us take any h G D(M, D*J\I) and any / G ©(M, DM). Suppose 
two different sets of fundamental solutions, say and 5^, exist. Starting from 

0={Ih-Ih,f) = {iS^-S^)Dh,f), 

uniqueness of the right fundamental solutions follows from the non-degeneracy of ( , ) , while 
uniqueness in the sense of left fundamental solutions can be deduced in a similar way. □ 

In analogy with the scalar case, we shall from now on call S = — S*" the causal propa- 
gator for the Dirac operator D and 5* = — S~ the causal propagator for D'. 

To conclude the section, we wish to underline that, up to this point, we have basically 
considered the Dirac spinor and cospinor fields as completely distinct objects. As it is usually 
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done in Minkowski space, however, we can define a well-behaved Dirac conjugation map mapping 
spinors into cospinors and vice versa. Furthermore, this conjugation turns out to map any 
dynamically allowed configurations into another one. 

Definition 2.10. We call Dirac conjugation matrix the unique matrix (3 G SL(4,C) such 
that 

P* = I3, i: = -haP-^ Va = 0,...,3, 

and furthermore il3n"'"fa > 0, being n timelike and future-directed. 

Starting from this object, we can define Dirac conjugation maps.- 

•t : E{DM) — > e.{D*M), = /*/?, 

: £{D*M) — > £{DM), h^ = (i'^h\ 
where * denotes the adjoint with respect to the Hermitian inner product on C^. 

Remarks. /3 is only unique, once a representation of the Dirac-Clifford algebra has been chosen. 
A direct inspection of the above identities shows that, with the definition of 7-matrices as in 
1^, P = -z7o and thus P = 

It furthermore follows in general that applying the Dirac conjugation twice gives the identity 
and that, as already anticipated, preserves the Dirac equations in the sense that D'f^ = (Df)^, 
Dh^f = (D'h)^ for any / G £{DM), h G £{D*M). 

3 Dirac fields: quantum point of view 

The aim of this section is twofold: on the one hand, we shall discuss the already available 
formulation of a quantum theory for Dirac fields on a curved background while, on the other 
hand, we shall show for the first time that the notion of Wick polynomials can be coherently 
introduced also in this scenario, giving rise to an enlarged algebra of observables. Particularly, 
we shall point out how all these topics fit into the framework of the locally covariant formulation 
of quantum field theory. 

To achieve our goal, we shall refer to an earlier work due to Araki |Ar70| . though we shall 
also profit from |Dim821 [Ko95l IHol991 IKranOl ISaVeOll IVerOU lE^Ml [Sa08] . 

3.1 The local algebras of fields and observables 

Although the first paper formulating a quantum theory of Dirac fields on curved spacetimes in 
the algebraic framework is |Dim82| , its underlying approach is slightly different from the one we 
shall employ, albeit it is ultimately fully equivalent. Particularly, in the aforementioned paper, 
the quantisation scheme calls for the choice of both a Cauchy surface S and initial data on S as 
building blocks of the quantum theory. We shall not dwell into the details of this method since 
we reckon that, in the spirit of local covariance, it is not best suited for our later purposes. 
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Although the standard paradigm in particle physics calls for the treatment of particles and 
antiparticles as distinct, albeit related objects, in this paper we shall, as it has been done by 
most of the authors mentioned in the introduction of this section, bear in mind the lessons from 
|Ar70j and, thus, we shall consider spinors and cospinors as part of a single entity, since it will 
turn out to be more convenient for our later purposes. 

On a practical ground, the building blocks of our discussion will be three. The first one 
arises out of the direct sum DM D*M, namely, 2) = 'D(M) = T){DM D*M), the space 
of compactly supported smooth sections of DM (B D*M with the usual topology, i.e., the one 
induced by the family of seminorms 

ll/IU = sup|/W(x)|, /gD, 

while the second is the map F : D — > 2) such that 

T{feh) = h^ep yfehev, (i6) 

being the Dirac conjugation introduced in definition 12. 101 Furthermore, in order to eventually 
impose the anticommutation relations, we need a third datum, namely, a sesquilinear form on 
= V {{DM D*M)^2^: let / = /a and /i = /ii /12 be two elements of V. Then 
(•,•) : 2)2 ^ C is defined as: 

{f,h) = -i{flSh,) + i{SMjh, (17) 

which is positive semidefinite, as one can infer with minor modifications either from lemma 
4.2.4 m [Sa08] or, with little more effort, from proposition 1.1 in [Dim82j : let us also note 
that (F/, F/i) = (/i, /). Furthermore, one can straightforwardly show that, using all the afore 
introduced tools, we can define the following algebra: 

Definition 3.1. We call algebra of fields the unital *— algebra 9'{M,g) generated by the 
identity and the abstract elements B{f) with f & D satisfying the following requirements: 

i) the map f ^ B{f) is linear, 

ii) B{Dfi D'f2) = for all /i /a E 

Hi) BiVf) = B{f)*, for all f (z D and with F defined as in (jl6p . 

iv) {B{fy,B{h)} = B{f)*B{h) + B{h)B{f)* = {f,h) , where the right hand side is given by 

m- 

In order to convince a potential reader that we are employing a sensible definition, let us 
first discuss some general properties. 

Remark. It is possible to recover the standard notion of spinor and cospinor quantum field 
starting from the i?-generators as follows: 
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the spinor arises as ip{h) = B{0 (B h). 



• the cospinor is given by ip^{f) = -B(/ © 0), 

Particularly, to be convinced of the self-consistency of such statement, one should notice that 
the spinor and cospinor fields are related due to property iii) and they respectively satisfy the 
Dirac and the dual Dirac equation of motion in the distributional sense thanks to ii). Finally, 



In order to better grasp the structure of 3'{M,g), one should realise that it is nothing but 
a topological *-algebra; this can be fully understood starting from an a priori different, albeit 
ultimately equivalent, perspective, namely, the so called Borchers-Uhlmann algebra for Dirac 
fields, which is explicitly discussed in |Sa08] . It is constructed as the following quotient 



where D" = D (^{DM © DAf *)^") denotes the compactly supported sections of the n-fold outer 
tensor product of DM © D*M and T)^ = C, whereas J is the closed *-ideal which arises out of 
the relations (ii) , (iii) , and (iv) in definition 13. 1[ It is hence generated by elements of the form 
Dfi © D'f2 with /i © /2 G D, by those of the form (/i © /a)* - r(/i © /a), and finally by those of 
the form / (g) /i + /i ©/ — (/, /i) with f,h £ D. 9"(M, g) can be equipped with a natural topology 
induced from that on This is tantamount to the request that a sequence fj = ®nfj,n 

is said to converge to / if and only if 

1. every fn in with respect to the topology of uniform convergence of all derivatives 
on a fixed compact set, 

2. it exists an E N such that fj^n vanishes for every n > N and for every j. 

In the forthcoming discussions it will sometimes also be possible to use a weaker version of 
3'{M, g) which is defined in the same way as 3^{M, g), but without including the Dirac equations 
in the construction of the ideal J; we shall refer to this case as the off-shell formalism. 

As a subsequent natural step, one can realise that the sesquilinear form (•, •) can be promoted 
to a genuine non-degenerate scalar product on the coset space D/ker 5 ©5"*, which, in turn, can 
be completed to a Hilbert space 'K with respect to the said scalar product. As a by-product, this 
entails the possibility to extend 9'{M,g) to a C*— algebra ^{M,g) representing the elements as 
bounded operators on itself. Following |Ar70] . we shall refer to this scenario as the assignment 
of the C*-algebra of Dirac fields ^{M,g) to the pair (J{, F). 

We have explained how to construct the local algebra in the case of Dirac fields; however, 
from a physical point of view, observables are required to commute at spacelike separations and 





(18) 
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the full 3'{M, g) does not fulfil such requirement. As a first step towards a definition of the 
algebra of local observables of a Dirac field, we can restrict our attention to 

3^even{M, g) = even subalgebra of 3'{M,g), 

which, e.g., can be defined as the subalgebra invariant under B{f) i— > —B{f) |Dim82t IHol991 
ISaOS]. The reason to choose such a subalgebra stems from the fact that any two elements of 
3^even{M^9) indeed commute for spacelike separations: 

Proposition 3.1. Let Ai, i € {1,2} he two elements of 3'even{M, g) which arise as finite 
linear combinations of smeared B{f) generators 

A = Y,B{f:,,i)---B{f:,^,,j 

n 

such that 

supp fnj and 

are spacelike separated. Then 

[Ai,A2] 



Proof. The proof descends out of two key observations: on the one hand we know the following 
relation between the commutator and anticommutator of four operators A, B, C and D: 

[AB, CD] = A{B, C}D - AC{B, D} - C{A, D}B + {A, C}DB . (19) 

On the other hand we know that, given two field algebra elements B{f) and B{g) with the 
support of / and g spacelike separated, condition iv) in definition 13. II together with the support 
properties of the causal propagator, proved in theorem 12.31 entail that 

BifyBih)+Bih)Bifr = 0. 

To conclude the proof, one needs to notice that, since only products of an even number of 
fields appear, the properties of the commutator allow to reduce [^1,^2] to a linear combina- 
tion of commutators, all of the form (I19p with AB and CD of the form B{f^ j_^)B{f^ j^) and 
-^(/n ji)-^(/n j^) ^respectively. This operation together with the requirement on the supports of 
the test sections defining Ai and A2 concludes the proof. □ 

With the restriction to 9'eveniM, g) we have been able to assure local commutativity. This 
criterion is, however, not sufficient to extract the observable elements out of 3"(M, g) and 
S'even{M, g) is thus Still too large to be a candidate for the algebra of local observables. To obtain 
such a good candidate, we have to take only so-called "gauge invariant" elements of 3^even{M, g) 
into account, c/., e.g., |Ar70t lDim821 [Hol991 ISa08| . and we denote the resulting subalgebra with 



U supp flj 
= . 
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A(M, g). Particularly, if we consider any A = Y,n ' " " ^(/n,2fc„) in S'eveniM, g), it is ly- 

ing in A{M,g) if and only if it is invariant under the action of any S € 5pino(3, 1); such an 
action is defined by a straightforward extension of the known one on DM and on D*M, first 
to DM © D*M and subsequently to arbitrary outer tensor products of the latter, such that we 
have a well-defined action on the test sections fn,i ® • • • fn,2k„ determining A. This definition 
of an observable is compatible with the product of algebra elements, and thus defines a subal- 
gebra of 3^even{M,g) in a well-defined way. In the rest of the paper we shall always work with 
3^even{M, g), though all our results can be applied to A{M,g) (and to the topological closures 
of the mentioned algebras) as well. 

It is remarkable that, in order to get to the definition of the various algebras introduced 
above, once a particular representation of the Clifford algebra has been chosen, the only other 
necessary datum is the geometry of the underlying manifold. This can be understood realising 
that, beside the Dirac bundles DM and D*M themselves, the overall analysis relies on the 
causal propagators S and S*, which are unique in a globally hyperbolic spacetime with spin 
structure. This apparently innocuous observation will play an important role in identifying the 
quantisation of the Dirac field as a particular locally covariant quantum field theory, as we will 
explain in the next subsection. 

3.2 Locality and general covariance 

In order to establish a connection between the previous discussion and the modern interpreta- 
tion of quantum field theory over curved backgrounds, it is mandatory to address the question 
whether the axioms of a locally covariant theory, as proposed by Brunetti, Fredenhagen, and 
Verch in |BFV03] are fulfilled for the above displayed algebraic quantisation of Dirac fields, and 
an affirmative answer has indeed been given in [Sa08j . We shall not dwell on a recapitulation of 
the precise definition of all the needed tools, e.g., the involved categorical notions here: instead, 
we choose to provide a short overview and we refer an interested reader to |BFVn3[ iBiTHM ISaM] 
for further details. That said, per direct inspection of the previous analysis, we can infer that 
the following axioms of a locally covariant theory are satisfied: 

1. It is possible to associate to every globally hyperbolic spacetime {M,g) with spin structure 
{SM,p) the corresponding *— algebra 3'{M,g) of fields in a unique way, once a global 
representation of the Clifford algebra has been chosen. 

2. To every map x which is an isometric and orientation preserving embedding of (Mi, gi) into 
{M2,g2) and at the same time maps {SMi,pi) to {SM2,P2) in a coherent and equivariant 
way (cf. Definition 2.3.1 of [Sa08]), one can associate an injective, unit-preserving *- 
homomorphism between the corresponding fields algebras 9'{Mi,gi) and 3'(M2, (72)- 

3. Let us choose two maps as above, namely, xi '■ (Mi, gi, SMi, pi) — > {M2, g2, SM2, P2) and 
X2 '■ iM2, g2, SM2, P2) gs, SM3, ps); then the following composition law is satisfied 
for the corresponding algebra morphisms 
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Let us note that the above axioms are also fulfilled in the off-shell formalism, i.e., for Dirac 
spinor fields not subject to the Dirac equations. We can, furthermore, add another two axioms 
in special cases: on the one hand, if we restrict the construction to 9'even{M, g), the axiom of 
Einstein causality is fulfilled on account of proposition I3.lt 

4. Consider two globally hyperbolic spacetimes with spin structure {Mi, gi, SMi, pi) and 
{M2, g2, SM2, P2) together with xi and X2 respectively, two embeddings into a third 
spacetime with spin structure (M^, g^, SM3, p^) of the aforementioned kind. Under the 
assumption that xi(-^i) and X2{M2) are spacelike separated in M3, it holds that for every 
Al e S'eveniMi,gi) and A2 G 3'even{M2, g2), [%i(Ai),ax2(^2)] = 0. 

At the same time, in the on-shell formalism - though not slavishly for the extended algebra 
we shall later introduce - the Time slice axiom (cf., proposition 4.2.22 of jSa08| ) holds: 

5. Let X '■ ~^ ^2 be a map between two globally hyperbolic spacetimes with spin structure 
with the properties already discussed. If a Cauchy surface of M2 is contained in x(Mi), 
then a-^ is an isomorphism. 

These axioms state properties of the full field algebras, but one can refine these statements 
and identify fields with a special behaviour under the maps x and a^, the so-called locally 
covariant fields pFVOSj IHoWaOl] . In fact, as discussed by Sanders in [SaOSj, the field B{-), 
and, thus, also the single fields ip{')^ V'^(')) are locally covariant fields. This entails that B{-) can 
be understood as family of continuous maps, indexed by spacetimes with spin structure iwH, 

5a/:D(M)^ J(M,<7), 

such that, given two spacetimes with spin structure Mi and M2 and a map x '■ ^1 ~^ -^2 with 
the properties discussed in axiom 1., one gets the same result by either building a quantum 
field out of a test section /mi € D(Mi) and then mapping this field to 3^{M2,g2) via or by 
mapping the test section to /a/j = X*(/a/i) G T)(M2) via the push- forward x* of X and then 
building BMiifAh) ont of it. On the level of maps, we thus have: 

° Bmi = ° X*- 

Similarly, one can identify certain observable composite fields as locally covariant quantum 
fields via specific choices of test section spaces, and, furthermore, the Wick polynomials we shall 
discuss later fit into the same framework as well. 

3.3 Spinors and Hadamard states 

The algebra A{M,g) C 3^even{M,g), is, at this point of our discussion, the best candidate to 
play the role of an algebra of observables for a free Dirac field theory. Unfortunately, this status 

■^We omit the other data determining a spacetime with spin structure in the remainder of this paragraph in 
favour of notational simplicity. 
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is far from being satisfactory because objects such as all the Dirac bispinors, the current in 
particular, and the (components of the) stress-energy tensor are not contained in A{M, g) or 
'Jeven{M-,9)- Since we want to consider these as genuine observables, the best option is to solve 
this problem along the same lines employed in the scalar case, namely, we shall suitably enlarge 
-^eveniM, g) to include all the wanted elements. Although reasonable and, to a certain extent 
natural, such idea comes with a price to pay, i.e., not all the well-behaved algebraic states for 
3'eveniM, g) are admissible for the extended algebra; in fact, we have to select only those with 
the suitable ultraviolet behaviour already possessed by the Minkowskian vacuum state. 

This is indeed not a novel problem and it has been tackled by several authors; the underlying 
philosophy is to characterise the admissible states imposing suitable physical conditions, such 
as finiteness of quantum fluctuations, and thus the possibility to employ these states to define a 
sensible expected stress-energy tensor [WaTTj . The translation of these ideas in a mathematical 
language leads to the notion of Hadamard states, which we shall now discuss in our frame- 
work. The available literature is immense and we point a reader interested in further details to 
|KaWa9H IE,ad96al IE,ad96b[ IKMK] for scalar fields or to |Hol991 IKraOni ISaVeOlj for a discussion 
related to spinors. 

As a first step and as main topic of this section, we shall proceed introducing the notion of 
Hadamard states for the whole 3'{M,g) and only later we shall restrict them to 9^even{^, d)- 
The already anticipated enlargement of the algebra to include all interesting observables of the 
free field will then be the core of a subsequent discussion. 

That said, henceforth, we shall consider a state a; to be a continuous, positive, and normed 
linear functional on 'J{M,g), such that 

u;{A*A) > VA G 3"(M, g), lo{1) = 1; 

since this algebra is generated, according to definition l3.lt by the abstract elements B{f) with 
f G T), every said state is uniquely determined by the set of its n-point functions, namely, 

u;n{fi,...Jn)=uj{B{h)...B{U)) 

where, due to the required properties a state, each ujn is a distribution on D". The bridge 
between the algebraic formulation of quantum field theory employed in this work and its usual 
Hilbert space description is in the non-trivial direction provided by the Gelfand-Naimark- Segal 
(GNS) construction (c.f., e.g., jHa92j ) which yields a representation of an algebraic state and 
a field algebra in terms of a Hilbert space state and operator valued distributions on the same 
Hilbert space respectively. Among all possible algebraic states, a distinguished role is played 
by the so-called quasi- free ones, whose n-point functions can be determined fully out of u)2. We 
shall focus on these, and, following |Ar70| . we recall: 

Definition 3.2. A state oj : 3^{M,g) ^ C is called quasi-free if, given any set of fi G D 

with i G {1, • • • , n}, uj {B{fi) ■ ■ ■ B{fn)) vanishes for odd n while for even n it holds 

n/2 

u{B{f,)---B{fr,))= J2 (-l)'""'n'^2(/.42.-i),/.„(2.))- 

TT„GS' 1 = 1 
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Here, S'^ denotes the set of ordered permutations of n elements, namely, the following two 
conditions are satisfied for Tin G S'^: 

iTn{2i - 1) < vr„(2i), 1 < i < n/2, 

TTn{2i - 1) < vr„(2z + 1). l<i < n/2 . 

Even though it is in principle possible to state the Hadamard property for general states 
|StrVeWo021 ISa08) . we will restrict our discussion to quasi- free ones and can thus state ev- 
erything on the level of two-point functions. In this context and on the level of single Dirac 
fields, two distinguished distributions appear: 

a;+(/,/i)=..(V(/i)V^(/)) and u:' {f ,h) = u (^^\f)^{h)) , (20) 

where / € D{DM) whereas h € 2)(D*M) and where both ip^{f) and V'(^) ^-re particular 
elements of 9'(M, g) as explained in subsection 13.11 Hence, it turns out that both uj^ and 
can be understood as distributions on T)[DM Kl D*M). 

We can now introduce the notion of Hadamard state and, as in the scalar case, it is remarkable 
and useful that this concept can be illuminated in two equivalent ways. The first one has recourse 
to the notion of wavefront sets |DuHo72| IHo90] . a concept which enables a refined formulation 
of a singularity structure of a distribution, and, to this avail, one should take into account 
that an a priori obstacle lies in the nature of the vector-valued distributions appearing in the 
context of Dirac fields. Particularly, since wavefront sets are more familiar in the context of 
scalar distributions, we need to specify how they can be defined for distributions with values in 
higher-dimensional spaces. To achieve this, it appears to be natural to define the wavefront set 
of a vector valued distribution as the union of the wavefront sets of the coefficients with respect 
to a (possibly local) basis-expansion and indeed this turns out to be an invariant concept due to 
the properties of scalar wavefront sets [Den821 IKraOOl ISaVeOlj . Specifically, we can define the 
wavefront sets of lo {x, y) = u) (x, y) E (x) E^'iy) as 

4 4 

^^('^^)= U U '^Vi^,!/)- 

A=l B'=l 

By defining wavefront sets in this way, we certainly loose information on the most singular "di- 
rections" of a vector-valued distribution. This information can be encoded in so-called polarised 
wavefront sets, as introduced in |Den82j and applied in [KraOOl IHol99| . Though of high math- 
ematical interest, such concept is of no use in our approach and we feel safe not to dwell into 
it since we would end up providing only shallow ideas. That said, let us state the first possible 
access to Hadamard states [Rad96al IKraOOl IHol99l ISaVeOT] : 
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Definition 3.3. A quasi-free state oj satisfies the microlocal spectral condition (^SC) 
and is thus called a Hadamard state if only if 

WF{U2) = [{x,y,k^,-ky) eT*M^^\{), \ {x,k,) ^ {y,ky), A;, o} . 

Here, {x,kx) ~ iu^ky) implies that it exists a null geodesic 7 connecting x to y such that k^ is 
coparallel and cotangent to 'j at x and ky is the parallel transport of kx from x to y along 7. 
Finally, kx>^ means that the covector kx is future- directed. 

Remarks. If a quasi-free state uj fulfils the /iSC, then lo^ possess the following wavefront sets: 

WF{uj^) = {{x,y,kx,-ky) eT*M^^\0, \ {x, kx) ^ {y, ky), kx^o] , 

where kx<0 states that kx is past-directed. 

An even stronger relation between the two distributions uj^ arises if we employ the anticom- 
mutation relation since it entails that 

u;+if,h)+ij-if,h) = -i{h,Sf) . 

By contrast, the distributions u}{ip{hi)tlj{h2)) and '^(V'^(/i)'0^(/2))) which, together with co^ 
determine uj2, have smooth integral kernels. For oj[ij){hi)^{h2)), this can be proved employing 
a symmetry argument already used in a similar way in |Rad96aj : due to the anticommutation 
relations, we have 

U{i,{hi)ij{h2)) = -0J (V(/i2)V'(/il)) ■ 

Hence, if {x,y, kx, ky) is an element of the wavefront set of the distribution on the right hand 
side of the previous equation, then (y, x, ky, kx) must lie in the wavefront set of the other one. At 
the same time, on account of the //SC, we know that W F {ijj{il:{x)il:{y))) is not invariant under 
the exchange of coordinates. This entails that WF {uj{il){x)'4){y))) = 0, hence, u;(V'(/ii)V'(^2))5 
and analogously u;('i/'^(/i)^^(/2)), possesses a smooth integral kernel. 

Although highly elegant from a mathematical point of view and thus very helpful in ab- 
stract proofs, the microlocal definition of a Hadamard state is neither the first one introduced 
chronologically nor the easiest one to cope with on the level of explicit calculations. In fact, as 
already promised at the beginning of this section, there is a different, more explicit definition 
of a Hadamard state via the so-called Hadamard form. For scalar fields, this has been rigor- 
ously introduced in |KaWa9l] . while for Dirac fields a similar concept has been proposed by 
|Ko951 rVer96j . To introduce it, we need the notion of a convex normal neighbourhood, which is 
an open subset of M such that any two points x,y £ can be connected by a unique geodesic. 
On any convex normal neighbourhood, we can introduce the smooth halved squared geodesic 
distance a{x,y), and, finally, formulate the following definition: 

Definition 3.4. A quasi-free state uj is said to be of the Hadamard form if and only if in 
any convex normal neighbourhood the distributions kernels of u>^ can be written as 

ij^ix, y) = ^±D'y {H^{x, y) + W{x, y)) , 
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where the index -y stresses that the dual Dirac operator D'y acts on the y-variahle, and the 
singular Hadamard distribution kernels can be specified as 

hHx, y) = + V{x. y) In (21) 

Here, U, V, as well as W are smooth bispinors and V and W can be expanded in powers of a, 
viz., 

oo oo 

Vix, y) = Y, VM^' VT' W{x, y) = Y^ Wn{x, y)a{x, y)", 

n=0 n=0 

where X is a reference length, and cr±e{x,y) = a{x,y) ± 2ie{T{x) — T{y)) + with e > 0. In 
the above formula, T is a time function, such that VT is timelike and future pointing on the full 
spacetime {M,g). 

We furthermore require H^ to be bisolutions of the spinorial Klein- Gordon equations up to 
smooth terms, i.e., 

P^H^{x,y) eE{DM®D*M), PyH^{x,y) e ^DM ^ D* M) (22) 
and demand that their difference is specified by the fundamental solution of P, viz., 

H+{f,g)-H'{f,g)=i{g,Ef), 
where f G V{DM) and g G D{D*M). 



Remarks. The existence of a time function T is guaranteed on any globally hyperbolic manifold 
|BeSa05llB^Sa06] . 

Furthermore, a completely satisfactory definition of the Hadamard form requires some more 
work to rule out spacelike singularities, to circumvent convergence problems of the series V 
and W , which are only asymptotic, and, finally, to assure that the definition does not depend 
neither on a special choice of the temporal function T nor on the employed convex normal 
neighbourhood. For further details and discussions of these aspects and the existence of states 
of the Hadamard form we refer an interested reader to jFNW8l| IKaWa9l| IKo95t IVer961 ISaVeOlj . 

To determine the so-called Hadamard coefficients U, V, and W, one has to exploit the equa- 
tions ([221). this point, we would like to stress a slight conceptual difference between Dirac 
spinors and scalar fields: in the case of scalar fields, the two-point function fulfils the Klein- 
Gordon equation in both entries, and this property is thus inherited by its singular Hadamard 
kernel up to smooth terms. Contrariwise, in case of Dirac spinors, (I22p does not follow straight- 
forwardly from the fact that the two-point functions lo^ fulfil the Dirac equations. In more 
specific terms, if we recall the definition of uj^ ()20p . we know that they fulfil 

D'^uj^{x,y) = DyUJ^{x,y) = 0. 
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Consequently, 

D'.D'y {H^{x,y) + W{x,y))=Q, 
DyD'y (F±(x, y) + W(x, y)) = -Py {H^{x, y) + W{x, y)) = 0, 

and, thus, both PyH^ and D'^D'yH^ are smooth. The smoothness of PxH^ does however, not 
follow automatically from these considerations, but has to be required or proven in a way similar 
to the one displayed in lemma 5.4. of [SaVeOlj . 

We shall explicitly discuss the computation of U, V and W in appendix IA.3[ To this avail, 
the following proposition will prove to be very helpful: 

Proposition 3.2. Let H^{x,y) be the Hadamard distribution kernels of a state introduced in 
definition\3^ Then {D^ - D'y)H^{x,y) = {Dy - D'^)H^{x,y) are smooth. 

Proof. The overall strategy calls for combining a deformation argument as devised in the ap- 
pendix C of |FNW8l] together with the so-called theorem of propagation of singularities (c/. 
theorem 6.1.1 in |DuHo72] ). 

That said, let us proceed in logical sequential steps and consider any Cauchy surface S ^ 
{M^g^u) of the spacetime we are interested in and let us choose an open neighbourhood of S, 
say Os, such that it is a causal normal neighbourhood of S, i.e., S is a Cauchy surface for Os 
and for each p,q (z Os such that p G J^{q), it exists a convex normal neighbourhood containing 
J~{p)r]J^{q). The existence of such sets in a globally hyperbolic spacetime and for any Cauchy 
surface E was first proved in |KaWa91j . 

The above mentioned deformation argument grants us that it is possible to construct an 
isometric, orientation and time-orientation-preserving embedding, say of Oe in causal normal 
neighbourhood Os' of a Cauchy surface E' of a second globally hyperbolic spacetime M' . Fur- 
thermore, one can engineer M' in such a way that, in the past of x(Os); it exists another Cauchy 
surface S" with a neighbourhood Os" which contains the image of a suitable neighbourhood of 
a Cauchy surface S'" in Minkowski spacetime under an isometric, orientation preserving, em- 
bedding X; and it is straightforward to extend x and x in such a way that they respect the spin 
structures. 

bmce on O2 X Os are constructed only out of the local geometric data via (I22p . it 
is possible to build a second pair which coincides with the push-forward under x of 
in x(Os) X x(Oe)- Furthermore, due to the propagation of the Hadamard form as proved in 
[FNW8H ISaVeOlj . H"^ are of the Hadamard form in 0^" x 0^" as well, and their pull-back to 
Os'" X Oe'" thus coincide with the Hadamard distribution kernels in Minkowski spacetime. 

Let us now consider 

= (D, - D'y)H^ 

and proceed to prove that these distributions have empty wavefront set. According to the 
above discussions, we can push-forward to M' and subsequently pull them back to the 
neighbourhood Os'" in Minkowski space in a well-defined way. In the Minkowskian region, the 
pushed-forward and then pulled-back versions of are identically vanishing and thus have 
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empty wavefront set, since the flat spacetime Hadamard kernel only depends on x — y due to 
translational invariance. To access the wavefront set of in the original spacetime region 
under interest, let us note that these distributions satisfy PxPyU^ = up to smooth terms and 
the same holds for their mentioned push-forwards and pull-backs, where the operator PxPy is 
properly supported, of real principal type, and homogeneous of degree 2 since it is the tensor 
product of two second order hyperbolic differential operators. From this it follows due to the 
contravariant transformation properties of wavefront sets under diffeomorphisms and the prop- 
agation of singularities theorem (see [SaVeOlj or chapter 8 of [Ho90j ) that the wavefront set of 
on Os X Os can only contain elements of the form 

{x,y,kx,0) or {x,y,0,ky). (23) 

Following a line of argument employed in the proof of theorem 5.8 in [SaVeOlj . we can infer that 
WF{u^) = in the following way: since are constructed as Dirac derivatives of and 
properly supported partial differential operators like D and D' do not increase the wavefront set, 
we know that WF{u^) C WF{H^). If we furthermore recall that specify the singular parts 
of uj^ and that these kernels have the " antisymmetric" wavefront set displayed in definition 13.31 
it follows that WF{u^) can not even contain elements of the form (|23p end are thus empty. □ 

As a result of the procedures described in appendix I A. 31 U and V turn out to depend only 
on the local geometry and the mass m, while the state dependence of io^ is encoded in W. This 
"universality" of the singularity structure of states of the Hadamard form allows for a locally 
covariant definition of normal ordering, as we will see in the next subsection. To this avail, it 
will be useful to compose the Hadamard distributions to a single object living on 2) (g) D, viz., 

Wi ® /2, /ii e h2) = {D'yH+) {hi, /2) - {D'yH-) (/i, h2) (24) 
= H+{hi,Df2)-H-{h,Dh2), 

where /i © /2, /ii /12 £ I?- Before we start working with Hadamard states, let us state the 
already anticipated and fruitful equivalence of the Hadamard form and the fiSC, which is a 
result due to [KraOOl IHol99[ ISaVeOT] : 

Theorem 3.1. Let us consider a state uj on 3'{M,g) with two-point function uj2- This 
satisfies the microlocal spectral condition, if and only if the distribution on T) 0D defined by 

f^h ^ uj2{f,h)-f){f,h), (25) 

has a smooth integral kernel, and, thus, uj^ are of Hadamard form. 

To conclude this section, we would like to mention a most notable property of Hadamard 
states: despite the well-known problem to fix a unique vacuum state for a quantum field theory 
on a generic spacetime, Hadamard states turn out to be "almost" unique since they are all locally 
quasi-equivalent [Ver94l IVer96l IDaHoOG] . This implies that locally the density matrix states on 
the Hilbert spaces obtained from Hadamard states via the GNS construction are all equal, and. 
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in more physical terms, that any two Hadamard states have a finite energy density with respect 
to one another. The latter statement is of course related to the expected stress-energy tensor, 
the main topic of the last section of this work. 

3.4 On the notion of Wick polynomials 

In the development of quantum field theory, a well-known obstruction arises whenever we con- 
sider the product of two fields, which, being distributions, cannot be safely multiplied unless 
special conditions are met. Since, as already anticipated, our ultimate goal is to enlarge the 
algebras under consideration to include observables such as the stress-energy tensor for Dirac 
fields, we are lead to tackle this problem. Like in the case of scalar fields, this results in the 
introduction of Wick polynomials and in the following we shall try to adapt an approach similar 
to the one discussed in the work of Brunetti, Duetsch, and Fredenhagen |BDF09| which in turn 
is related to further earlier works [BrFrOOl IHoWaOl] IHoWa02] . 

Unsurprisingly, in our scenario, there are differences to the above mentioned works due to 
the vectorial nature of our fields and their anticommutativity. This necessitates a treatment of 
Wick polynomials of Dirac fields on curved spacetimes on its own and we will thus develop them 
in this subsection since they have not been treated in the literature in the past. 

As already anticipated, upon enlargement of the field algebra to include Wick polynomials we 
have to restrict our state space to Hadamard states, which seems not to be a real loss since these 
are already distinguished and presumably the only "physical" ones for the algebras discussed in 
subsection 13.11 

As a starting point to define the extended algebra of fields, it will be more convenient not 
to start directly from 9"(M, g) or its subalgebras, though we shall consider the set 

oo 

n=0 

where the subscript A indicates that, for n > 0, one takes into account only antisymmetric 
elements, while D^^ = 2)° = C. Notice that it is required that a generic element F £ Q{M,g) is 
unambiguously determined by a finite sequence {F^")} of antisymmetric elements lying in D". 
This entails that, in a basis expansion with respect to £®(x) = Ea{x) E^{x), each element 
p{n) — -Fq"^..q^^£®^ • • • (X) £®" has antisymmetric coefficients, viz., 

©l,...,0fc,0fc4.l,...,0„ V^l' • • • ' "^^^ ^fc+li • • • ) Xn) 

= --^£l...,efc+i,0fe,...,0„(3^i'---'^fc+i'^fc,---,a;„) Vl<A;<n. (26) 

The set C(M, g) can be promoted to an algebra with respect to the following product which 
we shall henceforth indicate as -a'-, let F = {F^")} and G = {G^")} be two generic elements in 
e(M,5r), then 

{F-aG)^'''^= ^ a(f(p) , (27) 

p+q=n 
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where A is the operator of total antisymmetrisation such that F -j^ G is indeed an element of 
C{M,g). Specifically, A leaves invariant, while, for an arbitrary = /i o • • • (g) /„ with 
fi € D and n > 0, the antisymmetrisation reads 



A(/l ® • • • ® A) = ^ (-1)'""' ®---^fn 



where the sum is taken over all permutation^ 7r„ € Sn and A can be extended to T)^ by linearity 
and continuity. The algebra {Q{M,g), -a) can be interpreted as the algebra of functionals, in the 
sense of distributions on smooth sections, on the classical field configurations of Dirac spinors. 

The standard quantisation scheme is eventually realised changing the product -a into a 
suitable ^-product compatible with the anticommutation relations. The overall procedure, once 
a functional A : ^ C is selected, can be realized out of the map Fa : ^j^jQgg 
action on a generic element -F*-""^ of ID" (notice, here taken without antisymmetrisation), is 
required to be trivial if n < 2, whereas, for n > 2, 

n— 1 n „ „ 

TaFW =^ j df,{x,) J d^(x,)(-ir*+iAe»e^(x„x,)4';)..0jxi,...,x„) x 

i=i j=i+i^.j- 

X £®i(2;i) • • • £®^) • • • £®^) o • • • £®"(x„). (28) 

Here, A(xi,Xj) = A®*®^ (xj, x^) £0^ (xj) (g) £e^(a;j), with £e(a;) = E^{x) (B Eb{x), denotes 

the integral kernel of A, whereas the symbol £®*(xj) indicates that £®'(xj) must be omitted. 
On account of the regularity of C{M,g), we can safely define a ^-product as 

F *5 G = A as (F (g) G) , (29) 

where as is defined formal exponential! 



as = exp ( i-T^ 



here, Tg arises from (j28p if one inserts for A the functional 

3/10/2^51052 ^ S(/i 0/2,510 52) = -(/2,55i> + (5,/|,<7^), 
with and S^, being the causal propagators constructed in theorem 12.31 

Remarks. If we introduce a *-operation on (C(M, g),*s) via the straightforward tensorialisation 
of r (jl6p . the result is naturally isomorphic to the off-shell version of 3'{M,g) with its standard 



*Of course not all permutations employed are necessary, since in (|27[) the constituents will already be antisym- 
metric. The antisymmetrisation as defined here, however, is still valid and it constitutes the easiest way to write 
it without unnecessarily getting lost in combinatorics. 

^Since C{M,g) contains only finite sequences of test sections, the exponential series will always terminate after 
finitely many terms. 
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product. Particularly, B{f) G 3'{M,g) corresponds to / € Q{M,g) and the anticommutation 
relations on 9"(Af, g) correspond to 

f*s9 + 9*sf = iS{f,g) = {Tf,g), 

as follows by straightforward computation. The equations of motion can then be implemented 
by dividing out a suitable ideal. Since the Dirac equations will not be necessary in the following 
discussion, we will denote both the on-shell and off-shell algebras with (S(M, g),*s) and we shall 
consider them as being isomorphic respectively to 3^(M, g) and to its off-shell version. 

Up to now we focused on rather regular objects constructed out of D, but, alas, this does 
not suffice to reach our ultimate goals; as a matter of fact, we need to consider the spaces of 
compactly supported distributionqj as well: 

£'o = C, g''^ = £' (^(DM e D*M)^'^) , £' = £'^ 

The underlying leitmotiv is rooted in our interest in objects like Jj^^ dfi{x) f (x) -.tp^ {x)tlj{x) : 
which will in the subsequent discussion correspond to distributions like f{x)6{x,y), these are 
nothing but elements of £ " which are supported on the thin diagonal 

Diagn = {(xi, . . . , a;„) G M" j xi = • • • = x„}. 

Since this amounts to potentially ill-defined operations such as taking the product of distribu- 
tions at the same spacetime point, we cannot blindly extend (C(M, g),*^) (and equivalently 
9"(M, (7)) to incorporate these new objects into an enlarged *— algebra, but we have to require 
some suitable regularity conditions. 

Definition 3.5. We call extended set of functionals Q(.xtiM,g) the set containing finite 
linear combinations of compactly supported distributions F^") E £'" whose wave front set satisfies 
the following requirement 

WFiF'^''^) n (M" X (F^ uF")) = 

where y+ and V - are the closure of the future and the past light cone respectively in the fibre 
of the cotangent bundle at each point of M . 

In order to adopt this definition, we have to make sure that Cext(Af, g) can be made into an 
algebra and it is manifest that the product -kg is not up to the task since it would lead us to 
pointwise product of causal propagators, which is ill-defined due to their wavefront set. In order 
to avoid the aforementioned problem, we can replace -kg by which is nothing but (|29p with 

^Notice that elements of £' test compactly supported sections on D* M ® DM and not on DM ® D* M. This 
"dual" notation, as employed, e.g., in [Dim82] . is used to stress that D ^ £'. 
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S replaced by —2iS^ as defined in (p^ . This new ^-product is equivalent to the old one when 
this is well-defined, being 

F*HG = aH (a^ (^) *s «h (G)) , (30) 
OH = expTn, 

where Th is defined as in (j28|) upon inserting — i/2S for A. If we take into account the 
singular structure of Sj and the criterion to multiply distributions as devised by Hormander ( cf. , 
|DuHo72] or chapter 8 in |Ho90j ) . it turns out that the problem with -kg disappears since the 
pointwise product of the integral kernel of Sj with itself is a well defined distribution. 

Remark. The outcome of the preceding discussion is the introduction of {QextiM, g),-k}{). 
Recalling that {Q{M, g),-ks) has been isomorphic to 3'{M,g), we can reverse this viewpoint 
and just define the extended *-algebra 3^e.xt{M,g) = {Qexti^, g),*H)- Similarly, restricting the 
possible test sections and distributions taken into account, we can define the extended algebras 
3^even,ext{M, g) and Aext{M,g). That said, following slavishly the analysis of the scalar case in 
[HoWaOl] IHoWa02] . the product ([30]) rephrases Wick formula in the Dirac scenario. 

As already announced in the prior discussion, due to the form of the wavefront set of 
Hadamard states, we can extend them and only them to genuine states for 3'ext{M, g). Par- 
ticularly, this entails the standard paradigm according to which the product of two fields, say 
ip"^ {x)^p{y), should be regularised as 

■.^Hx)My) ■■ = i^\x)i;iy) + ^D'yH-{x,y) 

such that, for a Hadamard state lo, uj{:tl)\x)'ip{y) :) = (87r^)^^M^(x, y). At a level of expectation 
values, this can be equivalently seen as leaving Tp^{x)ilj{y) unchanged while io~{x,y) becomes 
uj~{x,y) + DyH~ {x,y). This somehow heuristic comment prompts the following: 

Definition 3.6. Consider a quasi-free Hadamard state lo, whose n-point function is indicated 
as ujn- One can define the regularised n-point function :ujn- as 

: : = = , if n is odd 

n/2 

:ujn- {xi,...,Xn) = X] 11(^2 -iD) (x^„(2j-i),2;^„(2i)) if u is cvcn, 

TTn&S'„ i = l 

with S) as in (p^ whereas the set of ordered permutations S'^ C Sn is the one introduced in 
definition \3.Sl 

Remark. As a straightforward consequence of the last definition, we can form expectation 
values of all elements in 3^extiM, g). Specifically, for any F = {F^""^} G S'e^tiM, g), 

n 
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is well defined due to the wavefront set properties of both the state and the allowed F. 

At this stage we need to point out that there are still some ambiguities in the employed 
definition of and thus in the definition of both and :uJn-', indeed, the reference length 
A necessary to construct according to definition 13.41 is in principle undetermined. This 
fact does, however, not hamper our analysis since different choices of A and thus of f) lead to 
isomorphic algebras. 

Lemma 3.1. Suppose we choose two different Sj, say Sji and Sj2, to construct the extended 
algebra {Cext{M, g),-ki{). Then the two resulting algebras {Cext{M, g),-ki{^) and {Qext{M, g),-*:^^) 
are isomorphic. 

Proof. Due to the properties of the Hadamard distributions H^, one knows that the difference 
d = S^2 — ^1 lias a smooth antisymmetric integral kernel. The two products and ■kH2 are 
related by a deformation. They are thus equivalent and the operator intertwining them can be 
realised as 

ad = exp (Trf) , 

where F is taken as in (|28p with d being inserted in place of A. Particularly, 

F-kH^G = ad {a'^i^F oQ^G) 

which is well-defined and holds true since a/fj ° cif^^ = ad and d has a smooth integral kernel. □ 

To finish the preparations for the final section of this work, we have to address a last issue. At 
the moment we are falling one step short from our ultimate goal since, to study the regularisation 
of the stress-energy tensor, one has to understand the treatment of differentiated fields. Hence, a 
small addendum to the above analysis is needed and we shall follow the procedure employed for 
scalar fields in |Mo03| . though adapted to our language. Thus, let us take a differential operator 
K on V{DM) of the form 

= ao + + ...Vf^, i2<cx), 

where 

and af ■■■^'=^5 for A; € {0, • • • , i?} are the coefficients of an element of 

Ti TM ®- - -®TM ®DM (g) D*M). 

k 

Notice that this class of differential operators encompasses both the Dirac operators and the 
spinorial Klein-Gordon operator which will appear in the expression of the stress-energy tensor 
and of its trace. In an analogous way we can choose differential operators K' on 1){D*M) 
and combine them with a K to operators K (B K' on D. If we now bear in mind definition 
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13.51 we realise that the extended set of fields is defined out of a condition on the wavefront set 
of its elements. Thus, in order to engineer any operator of the form K ® K' into the above 
discussion, we just need to recall a general result on wavefront sets (cf. Chapter 8 of |Ho90j 
or |DuHo72] ) according to which a partial differential operator which is properly supported, 
such that it maps compactly supported distributions to compactly supported ones, does not 
increase the wavefront set of a distribution it is applied on. Since the differential operators 
we are considering are properly supported, we can readily conclude that operators of the form 
K (B K' map Qext{M,g) to itself and that the previous discussion has already encompassed 
the treatment of differentiated fields. One could now prove several further properties of Wick 
polynomials of differentiated fields, but we will not indulge in this task since it will play no role 
in the forthcoming discussion and, furthermore, the results are by all means a straightforward 
extension, both as concepts and as technical proofs, of those discussed in [Mo03j for the scalar 
case. 

Before proceeding with the discussion of the stress-energy tensor, let us finally remark on how 
the extended algebra 3^extiM, g) fits into the framework discussed in subsection 13.21 Without 
going much into details we would like to point out that any coinciding point limits of smooth ob- 
jects constructed out of the Hadamard distributions are locally covariant, since the construction 
of depends only on the mass and the local curvature. As a result, all elements of 3'ext{M, g) 
which correspond to distributions with support on the thin diagonal are locally covariant fields. 



4 The stress-energy tensor of Dirac fields 

The aim of the section is to focus on the structure of the stress-energy tensor and to study its 
quantum properties. Particularly, we shall display that it is possible to introduce an improved 
tensor which is conserved also at a quantum level, although its trace acquires new and classically 
unexpected terms of geometric origin which lie at the heart of the so-called trace anomaly. 



4.1 The classical stress-energy tensor 

We start our analysis by revising the form and the properties of the stress-energy tensor for 
Dirac spinors in a classical framework. The Dirac equations (llip can be realised as the extremal 
of the unique action functional 



S 



d x^/\g\L 



d x\/\g\ 



M 



M 



A direct inspection of the above action shows us that, up to a total derivative term, it is identical 
to the more common expression 



S 



d^x\l\g\'4)^ (— 7^V;Ai + "T-V") 



M 
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We define the (Hilbert) stress-energy tensor by the usual procedure, i.e., 



- 2 5S 

An exphcit reahsation of this last identity in the case of spinor fields is much more involved 
due to the underlying orthogonal Lorentz frames whose explicit dependence on the metric must 
be accounted for. Nonetheless, a lengthy and, to a certain extent, tedious calculation, fully 
developed in |FoRo04| . yields 



T, 



'ixv = ^ (i^\[Pu)i^ - 'iphifii'-u)) - Lg^v, (31) 

where () denotes idempotent symmetrisation; one should also notice that, being the field free, 
the Lagrangian vanishes on shell. 

If we contract (j3ip with g^^^ , we end up with the classical trace 

T = g^'T^u = \ (V'f^.T^^ - V'W^V'r) - 4i = -mVV, 

where, in the second equality, we have evaluated the left hand side on shell by means of 
Hence, as expected, the trace vanishes on shell for conformally invariant, i.e., massless, Dirac 
fields. 

If we consider instead the covariant conservation, we need to calculate 



which vanishes on shell. 



4.2 The quantum stress-energy tensor: the problem 

In the next step, we would like to define the quantum version of the stress-energy tensor for 
Dirac fields. Since we have a well-defined notion of Wick polynomials at hand, it would be easy 
to just take the classical expression for the stress-energy tensor and replace the occurring field 
monomials with their normal ordered quantum counterparts. This way, one would easily get an 
element of AextiM,g) which, GNS-represented with respect to a Hadamard state, would be a 
well-defined operator valued smooth function. As we will shortly see, however, this procedure 
would not yield a meaningful object. To understand this, let us take a slight detour and think 
about the properties we would like a quantum stress-energy tensor to have. 

From the point of view quantum field theory over curved background, the most important 
entity to take into account as a guide in the search for a good quantum stress-energy tensor is 
of course the semi-classical Einstein's equation, viz., 

G^,{x)=SttGuj{:T^,{x):), (32) 
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where G^^, denotes the Einstein tensor R^^ — ^Rg^u, G is the gravitational constant and :T^u- 
is a suitable regularised expression for the quantum stress-energy tensor; this equation thus 
describes the back-reaction of the quantum field on the background. The legitimate question 
which now arises, is under which circumstances this equation makes sense at all. Regarding the 
form of the equation, we will restrict ourselves to point out that it can be derived by formally 
expanding a quantum metric and a quantum field about any classical vacuum solution of the 
Einstein's equation and computing the equation of motion for the expected metric while keeping 
only "tree-level" graviton contributions and "loop- level" quantum field contributions. Since one 
discards "loop- level" graviton contributions, the equation derived in this manner can only make 
sense as a model equation, or maybe for special states. We refer the interested reader to [FlWa96| 
and the references cited therein for an exhaustive treatment of this topic while we shall continue 
dwelling upon the properties of the quantities appearing in ([32]) . 

The first, and certainly obvious, observation is that we need a regularised expression of the 
stress-energy tensor to obtain a finite expectation value, i.e., a finite right hand side. The next 
observation is that the left hand side of ()32p is a classical and "sharp" quantity, while the right 
hand side is a probabilistic object. Such a situation can of course only be a sensible one if 
the fluctuations of the probabilistic quantity involved are small in comparison to the quantity 
itself, finite in particular. These considerations are exactly those which led to consider and 
select the states of Hadamard type as the physical and reasonable ones among the myriad of 
states available in quantum field theory on curved spacetimes. As a matter of fact, if we define 
normal ordering by means of the Hadamard singularity and we evaluate the normal-ordered 
(and smeared) stress-energy tensor on Hadamard states, we automatically get a quantity with 
finite fluctuations. This stems from the fact that powers of the Hadamard bidistribution are 
again well-defined bidistributions, thanks to the special form of the Hadamard wavefront set. 
Regarding quantitative statements about the fluctuations of the expected stress-energy tensor, 
it seems that in general a priori statements are not possible and one has to look at solutions of 
the semiclassical Einstein's equations to a posteriori compute the fluctuations on these solutions 
and inspect to what extent these are to be trusted. 

It will prove helpful to realise that uj{: T^u{x) :) for a : T^uix) : given as some linear combi- 
nation of the previously defined Wick monomials and evaluated on a Hadamard state to can be 
equivalently expressed as 



where D^j^ is some (bi) differential operator specified by the choice of linear combination of Wick 
monomials in the definition of :T^jj{x) :, Tr denotes the trace over spinor indices and we refer 
the reader to appendix |A] for the explanation of the possibly unfamiliar notations that arise 
in the context of bispinorial entities and will be extensively used in the following. As we have 
already remarked, the obvious choice of expression for the stress-energy tensor in terms of Wick 
monomials will not turn out to the best one. In terms of the above defined differential operator. 




(33) 



38 



this means that the canonical version derived from the classical stress-energy tensor pip . 

is not well suited for defining a sensible uj{:T^y{x) :). 

Since we have assured ourselves that a right hand side of (j32p obtained by expressing uj{: 
Tfj^uix) '■) as ()33p is in principle well defined, we could seek for additional physical and consistency 
requirements that lead to a potential refinement of that procedure, i.e., to a sensible choice of 
D^^. Pursuing a comparable aim, Wald |Wa77l IWaTSj has set up five axioms that a meaningful 
expected stress-energy tensor should fulfil. These proved to be a valuable tool in a posteriori 
legitimating known stress-energy tensor regularisation schemes in curved spacetimes and stating 
to which extent they may differ from one another without raising doubts about their validity. 
For the convenience of the reader, we list them in the following. 



Definition 4.1. We say that uj{:Ti^i^{x) :) fulfils (the strong version of) Wald's axioms if it 
has the following five properties: 

1. Given two quasifree states loi and L02, such that uj^{x,y) — 1^2^(2;, y) is a smooth bispinor, 
uJi{:T^uix):) - u;2{:T^uix):) is equal to Tr D'^^^'J' {lj^ {x,y) - uj^ {x,y)) 

2. uj{:T^y{x):) is locally covariant in the following sense: let 

X : (Mi,gi,5Mi,/Ji) ^ (M2, 52, pi), 

a^:3'{Mi,gi) ^3^{M2,g2) 

as in subsection \3.^ If two states uji and 002 on 3^[Mi,gi) and 9"(M2, 52) Q^re related by 
(jji = (JJ2 ° ot-^, then 

uJ2{:Tf,^u2{x2) ■■) = X* i^^ii-T^iuiixi):)) , 
where x* denotes the push-forward of x in the sense of covariant tensors. 

3. yt'uj{:T^^{x):)=0. 

4. On Minkowski spacetime and in the Minkowski vacuum state coMink, ^Mink{'-T^u{x) '■) = 
0. 

5. uj{:T^jj[x):) does not contain derivatives of the metric of order higher than 2. 

ijo{:T^y{x) :) is said to fulfil the reduced version of Wald's axioms, if only the first four statements 
hold. 

Wald has originally stated the axioms for scalar fields, while the version we give here is 
modified to be suitable for Dirac fields. We reckon that a few comments both on the origin 
and on the meaning of the single axioms might be helpful for a potential reader in order to 
understand their relevance: 
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1. In a given Fock-representation of the quantum field, the non-diagonal matrix elements of 
the formal unrenormalised stress-energy tensor operator in the "mode basis" are already finite, 
because their calculation only involves "finite mode sums" , while the calculation of the diagonal 
matrix elements involves "infinite mode sums" (WaTTt IWa78j . To regularise the formal stress- 
energy tensor operator, it is therefore only necessary to subtract an infinite part proportional 
to the identity operator, thus leaving the non-diagonal matrix elements unchanged. Axiom 1 
amounts to require such a "minimal" regularisation. This axiom is also related to so-called rel- 
ative Cauchy evolution of a locally covariant field [BFV031 ISa08| ; since the functional derivative 
of the relative Cauchy evolution involves the commutator with the stress-energy tensor oper- 
ator, one could reformulate this axiom on the operator level requiring that any regularisation 
prescription yields the same relative Cauchy evolution. If we consider Hadamard states, such 
that uj^{x,y) is locally given by —Dy{H~{x,y) + Wj(x,y))/(87r^) for i = 1,2, the requirement 
is equivalent to demanding that the differential operator used in ([33l) is given by plus a 
term which does not influence the state dependence of uj{:T^y{x) :). 

2. Taking the locality principle of quantum field theory and the covariance principle of general 
relativity seriously, we would like to have a a;(: T^y{x) :) which describes the back-reaction of 
the quantum field on the spacetime in a local and covariant way. In fact, this axiom seems to 
have been an inspiration towards the formulation of locally covariant QFT, as described in the 
seminal paper |BFV03| . 

3. This axiom basically points out a necessary condition for the well-posedeness of the semiclas- 
sical Einstein's equations; namely, since the geometric left hand side of ()32p is conserved due to 
the Bianchi identities, also the right hand side should vanish under the action of the covariant 
derivative. 

4-- It is a sensible prerequisite of any regularisation scheme for a field theory on a curved 
background that it should be possible to read it as an "extension" of the standard normal 
ordering in Minkowski spacetime, but there are good reasons to skip this axiom, c/., |Fu89] and 
the second remark after theorem 14. 1[ 

5. Wald originally proposed this axiom in a rather technical and more strict way |Wa77t IWa78] , 
essentially requiring that Lo{:Tfj_y{x) :) does not depend on derivatives of the metric of order higher 
than the first. The underlying motivation is rooted, on the one hand, in the request of well- 
posedness of the Cauchy problem for the Einstein's equations even with a non vanishing source 
and, on the other hand, in the need for a sensible "classical" limit of the semi-classical Einstein's 
equations (see the enlightening discussion in Wald's original paper |Wa77l rWa78j ) . Wald himself 
had realised, however, that the strict version of this axiom could not even be satisfied in the 
classical theory and has, thus, proposed the weaker one stated here. Unfortunately, further 
examinations have revealed that even this weaker version does not seem possible to fulfil in 
massless theories without introducing an artificial length scale into the theory; therefore, the 
axiom has been discarded. We still believe, however, that it could be fulfilled, though only under 
special circumstances. We shall comment on this issue at a later stage of the paper. 
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Using these axioms, Wald could prove that a uniqueness result for a;(: T^y{x) :) can be 
obtained. The first two axioms already imply that the results from two different sensible regu- 
larisation schemes can only differ by a local curvature tensor. The third and fourth axiom then 
imply that this local curvature tensor is conserved and vanishes if the spacetime is locally flat. 
Requiring that this term has the correct dimension of m^, the possible tensors are presumably 
only the ones obtained by varying a Lagrangian of the form 

\ \m? J \ J \ J J 

with respect to the metric, with some dimensionless functions Fi{x). Requiring suitable analyt- 
icity properties with respect to the curvature tensors and m, in |HoWa05] . it has been shown 
that the only possibilities are m^g^i, {Fi = ij^, m?G^y {Fi = x) and the three local curvature 
tensors I^,^ {Fi = x^), J^i, {F2 = x), K^y (F3 = x), cf. appendix lAl for their full expression. In 
fact, we will later show that changing the scale A in the regularising Hadamard bidistribution 
amounts to changing lo{: T^y{x) :) exactly by a tensor of this form and, furthermore, the 
attempt to regularise Einstein-Hilbert quantum gravity at one loop order automatically yields 
a renormalisation freedom in form of such a tensor as well |tHoVe74] . Having in mind how 
the semiclassical Einstein's equation may be derived, these two arguments are of course related 
by means of internal consistencylf] Using the Gauss-Bonnet-Chern theorem in four dimensions, 
which states that 

j df,{x)R^yprR'''''" - 4.Rf,,R^"' + R^ 

M 

is a topological invariant and, therefore, has a vanishing functional derivative with respect to 
the metric |A1951 ltHoVe74] . one can restrict the freedom even further by removing Kfj^i^ from 
the list of allowed local curvature tensors. 

4.3 The quantum stress-energy tensor: the solution and its trace anomaly 

We now seek to exploit the above axioms in order to specify a sensible choice of differential 
operator D^^^. Looking at our proposed regularisation procedure (133]) . the first obstacle to over- 
come seems to be the covariant conservation axiom. As we have seen above, conservation of the 
stress-energy tensor in the classical case is a direct consequence of the equation of motion. Since 
we are regularising by subtracting from the two-point function the Hadamard bidistribution, 

term proportional to the metric is not allowed if one seeks to fulfil the third axiom. As we will see later, 
however, it does not seem to be possible to fix this term in a way that is compatible with analyticity in m. 
Furthermore, the results of Hollands and Wald regarding the restriction of the possible regularisation freedom 
by demanding analytic dependence on curvature and mass have only been obtained for scalar fields. Since the 
stress-energy tensor for Dirac fields is an observable and thus still a "scalar" field, their results can be, nonetheless, 
presumably extended to this case. 

*In fact, at least in the case of scalar fields, the combination of the local curvature tensors appearing as the finite 
renormalisation freedom in |tHoVe74] is, up to term which seems to be an artefact the dimensional regularisation 
employed in that paper, the same that one gets via changing the scale in the regularising Hadamard bidistribution. 
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which is in general not a solution of the Dirac equation(s), -D^"" applied to the thereby obtained 
smooth bispinor will in general not yield a conserved quantity. A viable solution to such prob- 
lems calls for the modification of the classical stress-energy tensor (j3ip by terms which vanish on 
shell, while, at the same time, they help restoring covariant conservation on the quantum side. 
As in the case of scalar fields |Mo03j . it seems that the only possible option is to add multiples 
of the Lagrangian to the classical expression of the stress-energy tensor. We thus propose the 
following classical stress-energy tensor as a starting point 



and we look for a c € M that yields a sensible uj{:T^y{x) :). This is tantamount to the choice of 
the following differential operator for the point-splitting process (|33p 



Before proceeding to prove that there indeed exists a suitable choice of c, we would like 
to anticipate another result, which can be easily understood from the aforementioned line of 
argument. Let us remember that there is another property of the classical stress-energy tensor 
stemming from the equations of motion: it has vanishing trace in the massless (and therefore 
conformally invariant) case. Following the above discussion, on the one hand, it seems that one 
might need to give up this property at a quantum level, while, on the other hand one, one could 
still hope that the choice of c also provides a vanishing trace. Alas, it will turn out that this is 
not the case. One can only fix c in a way such that ijo{:T^jj{x) :) has vanishing trace for m = 0, 
but conservation is inevitably spoilt. Since we have already realised that conservation is indeed 
an essential requirement for the right hand side of the semiclassical Einstein's equations, we will 
have to accept that g^"Lo{:T^j^{x) :) is not vanishing in the massless case. This goes under the 
name of trace anomaly. 

Theorem 4.1. Let Am = 2exp(| — 27)771"^ for m ^ and Am arbitrary for m = 0, where 7 
denotes the Euler-Mascheroni constant, choose the Hadamard bidistribution to be the one with 
A = Am and let uj{:T^,y{x) :) be defined as in (33\). with the differential operator D^jy = D^^^ 
defined as in [34\ )- Then uj{:T^y[x):) fulfils the reduced version of Wald's axioms. Furthermore, 





(34) 
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it exhibits the following trace (anomaly) 

r^i-.W :) = - ^ (^i^^ + i^UR - ^^R,.R^^ - ^,R..,rR^^-) (35) 



^ ^ Ic^uprC""'^ + 11 (r^.r'^'' -Ir']- qbr 



28807r2 ^^2^^^^^"^ ' V 3 
1 f m'^ m?R\ ^ v^i^^r, NT 



= Tr 



Proof. We begin by computing V^uj{:T^uix) ■) and g^^ ijj{:T^jy{x) :), leaving c unspecified for the 
moment. Applying Synge's rule and taking into account that [C^jy,7^] = and [f/J^.^] = (c/., 
appendix |A|) , we get 

Sn''V^Lo{:T^,{x) = V^Tr y)] = Tr [(V^ + 5^ V^^^.lx, y) 

\{ {Ou'^w - V,) (D^D; + P,) + ^7.^;(Py - Px) 

-| (5:;' V,, + V,) {D'.D'y - Py) } W{x, y) . 

Remembering that —D'y{H~ + W) is the local two-point distribution of a state, it follows that 
H^+W is subject to the distributional differential equations D'^D'y{H~+W) = = Py{H^+W). 
Thus, we can safely replace W in the above equation by —H~, since every appearing term 
involves one of the two aforementioned differential operators. Such a procedure yields 



87r2v^u;(:T^^(x):) =Tr 



(v. - <7:;'V,,) {D'.D'y + Py) + + m) (P.. - Py) 

<'V,, + V,) {Py - D',D'y)]H~{x,y) . 



Now we can insert the various coincidence point limits of the differentiated Hadamard bidistri- 
bution computed in proposition lA.il of appendix lAl to obtain 

87r2v^L^(:r,,,(x):) = -{I + Qc)Tr[Vi{x,y)],,. 

For the trace we use both the insights on the parallel transport of gamma matrices from appendix 
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1X1 and the arguments already employed in the computation of the conservation to get 



^7:''g^'-u{:T,,{x) :) = g^^Tr [Z)^,(x, y)W{x, y)] 



Tr 
Tr 



2c + {D'.D'y - Py) + ml?; I W{x, y) 



2c + 



1 



2 , ,D',D'y - Py) H- {x,y)+ mD'yW{x, y) 

= -Q{Ac + l)Tr[Vi{x,y)]+mTr [D'yW{x,y)\ . 

If we look at the above two results and we use the data from appendix [XJ we realise that, 
since Tr\Vi{x^y)] is in general neither vanishing nor constant, we need to set c = —1/6 to 
assure conservation, thus yielding the asserted trace and, particularly, the trace anomaly in the 
massless case. 

Let us now proceed to check the validity of the first two axioms. From the above two 
calculations we can extract that the term we have added to the canonical differential operator 
Z?^"" to achieve the looked-for differential operator L)^^ contributes to oj{: T^jy{x) :) a term 
proportional to 

g^,,Tr [{D'^D'y - Py) H-{x,y)] = -l2g^,Tr[Vi{x,y)]. 

This term is clearly independent of the state, i.e., thus, axiom 1 holds for our regu- 

larisation scheme and two Hadamard states oji,uj2- Furthermore, since the Wick monomials are 
locally covariant quantum fields as discussed in the previous section, the same holds true for 
:T^y{x) : and axiom 2 is straightforwardly fulfilled. 

What now remains to be shown is the vanishing of lo{:T^i,{x) :) on Minkowski spacetime and 
in the Minkowski vacuum state, provided we choose a suitable scale A in the definition of H^. 
In this setting, we have uj~^j-^j^{x,y) = — y)/4, with the scalar two-point function 

^2 Min^ |Scha95j . The latter can be specified as 



UJ. 



2, Mink 



{x,y) = lim 



4m 



0+ (47r)V2cTe(x,y) 



Ki{myj2a^{x,y)), 



where Ki denotes a modified Bessel function, the expression is to be understood in the sense of 
analytic continuation for negative values of the geodesic distance, and the limit specifies how to 
approach the branch cut of the squared root |Mo03j . Expanding this in terms of a (suppressing 
the e for simplicity) we have 



2,Mink 



;x, = ^ - + - + -a + Ha)a^ In ^ ) + 



m 



1 + 



'Recall that, according to our definition, uj is of "positive frequency type". 
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where the /j appearing in this paragraph are smooth functions. In Minkowski spacetime, the 
Dirac Hadamard bidistributions are simply the scalar ones times the unit matrix. This is 
not surprising since, as visible in appendix[Al the nontrivial matrix part of the Dirac Hadamard 
coefficients stems from the curvature of the spin connection, which vanishes in flat spacetimes. 
We thus have 

8-^^M.n.(-> = + + + 1- (J) } 

If we take into account this singular part, a short computation yields 

8n^WMink{x, y) = \-Y (in ( ^ J " ^ ) + ( ^ ( 2 ) " ~W) " ^ •^'^'"^ j ^ 

and we can straightforwardly compute 

, ^ , , , / /e^^m^A^X 7 
u)M^nk{-T^,u[x):) = — —g^i, (In ( ^ I - - 

which vanishes for A = Am = 2exp(| — 2j)m~^. In the massless case, both H^^^^ and uj2^Mink 
only consist of the cr~^ term, such that iOMink{-T^u{x) :) is trivially vanishing, independent of 
any scale A. □ 

It is interesting to notice that the above analysis seems to suggest that the common idea 
which associates the emergence of the trace anomaly to the conservation of the stress-energy 
tensor is somehow inappropriate since such anomaly seems rooted in the loss of the equations 
of motion at a quantum level alone. As a matter of fact the needed modification of the classical 
stress-energy tensor does not evoke the trace anomaly, but it only modifies it. 

To conclude the section, a few further remarks are in due course: 

Remark. Our result for the trace anomaly coincides with previous ones obtained by means 
of gravitational index theorems |ChDu79] and point-splitting techniques [Chr78j . Both ap- 
proaches have made use of the DeWitt-Schwinger expansion, which can not be defined rigorously 
|Wa79t [BrQt86l [Fu89l IMo03] . Moreover, even if this expansion reproduces the Hadamard singu- 
larity structure, it seems that calculations are much shorter if one expresses it directly through 
the Hadamard series. The previous attempts had, however, one advantage, namely, they had 
expressed the expected stress-tensor as the functional derivative of a (diffeomorphism-invariant) 
effective action, such that the result has been manifestly conserved. The nice discussion which 
follows definition 2.1 of |Mo03j gives an explanation of how the extra term in our derivative 
operator D^^, can be understood in this context. 



Remark. By changing the scale A in the Hadamard bidistribution H~ to A', one modifies 
the (definition of the) smooth part W by 2V In X' / X, such that a;(:T^;y(x) :) changes by a term 
proportional to Tr[D'^^,V]. Since we know from proposition 13.21 that D'yV fulfils both Dirac 
equations of motion, we can a priori deduce that this term is automatically conserved and 
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furthermore traceless in the conformal case. Thus, it follows that both the determination of the 
correct c to be inserted in D^^^ and the trace anomaly are independent of the scale A. Even if we 
already know the properties of Tr[D^^V] beforehand, it is enlightening to calculate its explicit 
form. The result is 

4 2 1 

TTl TTl 1 

TAO'l^uV] = ^5/.. - -^G^, + - (V - 3V) , 

where the linear combination of 1^^ and J^y appearing in the above formula is traceless, c/., 
appendix VK\ This term is well within and even exhausts the regularisation freedom discussed 
after definition 14.11 

Furthermore, one should recall that, in the massive case, A had to be fixed in terms of inverse 
powers of m to assure vanishing of oj{:T^y{x) :) in Minkowski spacetime. Hence, if one demands 
continuous dependence of u}[:T^y[x) :) on m, it does not seem possible to fulfil the third axiom 
of definition 14.11 in this way. 

Remark. Even if we are able to fix the scale A, we could in principle still add multiples of 
rin?G^y, I^y, and J^y to a;(: Tfj,y{x) :) without spoiling the validity of the first four of Wald's 
axioms, though possibly modifying the OR term in the trace anomaly. This freedom can be 
derived in a more general sense, by already viewing all Wick products as being uniquely defined 
only up to terms depending on the mass and the local curvature, where the possible regularisation 
freedom is partially restricted by suitable consistency conditions, e.g. "Leibnitz rules". This 
approach has been developed and pursued successfully by Hollands and Wald in |HoWa05| and 
has the advantage to also encompass interacting fields. A treatment along their lines will, as 
already remarked in subsection 14.21 presumably yield the same renormalisation freedom for the 
stress-energy tensor like the one found here, as it happens in the scalar case [ HoWaOS] . 

Furthermore, the arising of the conserved local tensors I^y and J^y puts us in the position 
to understand why the fifth of Wald's axioms in definition 14.11 is problematic. Let us consider 
the massive case, where we can fix the scale A. Since Ifj_y and Jfj_y contain terms involving 
fourth order derivatives of the metric, we may hope to cancel terms of that type occurring in 
uj{: T^y{x) :) by adding a fixed linear combination of those two tensors. In the massless case, 
however, there is, up to our knowledge, no physically sensible way to fix A. Therefore, one has 
no control on the multiples of I^y and J^y occurring in Lo{:T^y{x) :) and thus no way to cancel 
them. 

Nonetheless, there are scenarios where the situation with respect to the mentioned axiom 
is not that pernicious. On cosmological, i.e., Friedmann-Robertson- Walker backgrounds, the 
semi-classical Einstein's equations ()32p can be reduced to an equation for the traces of both 
sides plus a conservation equation for the right hand side [DFPOSj . Thus, it seems that one has 
the chance to fulfil the fifth axiom for both massive and massless fields in this simplified setting, 
since a change of scale does not add fourth order derivative terms to the trace of the expected 
stress-energy tensor. In fact, as already explained in the introduction, this observation has been 
used in [DFPOSj to obtain stable solutions of the semi-classical Einstein's equation at late times. 
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5 Conclusions and outlook 



We have extensively discussed the structure of free Dirac fields both at a classical and at a 
quantum level. While, in the first case, we have mostly reviewed standard approaches, in the 
latter scenario we have achieved a twofold goal. Particularly, we have started the discussion of 
quantised Dirac spinors by exploiting the selfdual framework introduced by Araki which treats 
spinor and cospinor fields as a combined single object and allows to formulate the quantisation 
procedure in a locally covariant way. This step has been fully undertaken by Sanders for the 
first time and we have recalled the essential steps and features of this construction. Employing 
already available and known properties of Hadamard states, we have subsequently been able 
to introduce the extended algebra of Wick polynomials, the topic of section 13.41 and the first 
of our main results. As a second one, we have shown that, as in the scalar case, a physically 
sensible definition of the stress-energy tensor for Dirac fields on a curved background in terms 
of Wick polynomials is indeed possible with just one caveat: one has to add to the classical 
expression a suitable term which vanishes on-shell and hence does not alter classical dynamics 
to obtain a conserved stress-energy tensor on the quantum side. Some new insights on Diracian 
Hadamard forms have constituted a prerequisite of this result, while one of its consequences is 
the emergence of a non- vanishing quantum trace of the stress-energy tensor, even if its value 
at a classical level is zero in the conformally invariant case. This result, which goes under the 
name of trace anomaly, has been previously known, but only as a result of formal calculations; 
it is thus derived here rigorously for the first time. 

On the overall, we reckon that this paper accomplishes also a further task, namely, it adds 
the insight that it is possible, interesting, but by no means straightforward to recast many of 
the already known rigorous results for scalar fields also for the spinor ones. Furthermore, our 
analysis opens several interesting questions to be tackled in future lines of research: the first 
one, which arises out of section 13.21 concerns the possibility to prove the time slice axiom for 
the extended algebra of fields (as well as for interacting field theories) in the scenario considered 
in the paper. If one follows the path paved in the scalar case in |ChFr08j . a positive answer 
seems definitively within our grasp. A further interesting problem originates from section 13.31 
in which Hadamard states are introduced and discussed; the Hadamard coefficients appearing 
in the singularity structure of such states are smooth bispinors and the question arises if their 
most remarkable feature in the scalar case, namely, their symmetry as proved in [Mo99t iMoOOj . 
also appears in the spinor scenario. Such a property would be desirable since, for example, it 
would lead to many simplifications in the demanding calculations necessary in the construction 
of the conserved stress-energy tensor. Although there are hints pointing towards this direction 
(see also |SaVe01j ) , we are far from a complete proof of such a symmetry and we thus feel this 
would be another rather interesting problem to tackle in the very next future. 

Besides these rather formal lines of research, our results have also some remarkable conse- 
quences at a physical level. On the one hand we are now ready to answer the question posed in 
the introduction on the robustness of the results in |DFP08| : preliminary considerations seem 
to point towards this direction, though we leave a definitive answer to a future analysis. On 
the other hand, since our approach allows us to control the behaviour of free Dirac fields at 
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a cosmological level, it is interesting to point out that free or perturbatively self-interacting 
fields with half-integer spin in cosmology arise in many models, such as baryogenesis through 
leptogenesis, where they often play a pivotal role. In these scenarios there are still many open 
questions to be answered and it seems that, often, the role of spacetime curvature effects are 
a priori discarded as negligible. Our experience suggests that this approximation might be too 
crude and, therefore, we would like to investigate these models in more detail in the framework 
of quantum field theory in curved spacetimes with the hope that such an analysis might lead to 
new and interesting physical consequences. 
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A Useful tools and necessary calculations 

The aim of the appendix is to recollect, to clarify, and, occasionally, to also prove useful formulas 
which are needed in the main body of the paper and which are subject to potential ambiguities. 
These are often perniciously leading to potentially grievous sign mistakes or misunderstandings 
of a sort which we wish to hold off from a potential reader. 

A.l Notations, conventions, identities 

As a starting point, we would like to recollect our basic conventions regarding some symbols, 
whose exact definition often varies among the literature. In accord with section [21 we work 
with spacetimes thought as four-dimensional, Hausdorff, smooth manifold^^, endowed with a 
Lorentzian metric g^i, with signature (— , -|-, +, -|-). At the same time, other notable geometric 
quantities, namely, the Riemann and the Ricci tensor as well as the Ricci scalar, are defined via 
their components as follows 

where Va are the components of an arbitrary covector; the extension to vectors and tensors of 
higher rank is then straightforward. As a last remark, we underline that the Riemann tensor 
possesses the symmetries 

Rafi-yS — Rfia'^S — RajBS-y — R'ySa/B 

^"As proven in [Ge68) . second countability is automatically fulfilled for a four-dimensional, Hausdorff, smooth 
manifold of Lorentzian signature. 



48 



and fulfils 

Finally, we define the Weyl tensor as it is usually done by the following expression 

1 1 

CaPyS = Ra/S-yS — g {9aSgf3^ + ga^Qps) R~ 2 (d/SsRa-y + g/S-yRaS + gaSRpy + ga-yRfSs) ■ 

Similarly, we need to cope with geometric quantities related to the spin structure, introduced 
in definition 12. 4i Most of these are constructed out of the so-called 7-matrices which satisfy the 
standard anticommutation relations (jH), i.e., {7/^,71/} = "^g^u- Our choice of the metric signature 
entails that the 7-matrices are different from the standard ones employed in quantum field theory 
books by an overall multiplicative factor iti. Consistently also with ([2]), we stick to +i and, 
therefore, the Dirac operator appearing in the Dirac equation for spinors becomes D = — y + m, 
whereas the operator nullifying a dynamically allowed cospinor is D' = Y + m. 

That said, apologising in advance for assigning the letter C to two different objects, we define 
the components of the curvature tensor C of the spin connection as 

where Va are the components of an arbitrary cospinor; as previously, the extension of this 
definition to spinors, also of higher rank, as well as of that for the Riemann and the spin 
curvature tensor in presence of mixed spinor-tensors, is straightforward. It follows from lemma 
12.11 that the relation between the two curvature tensors is 

^ B^lu — -^^tivapl c7 B- 

Thus, C possesses the symmetries 

CaB/iu = ~CBAfiu = ~Cabu^- 

We also use the notational convention that a matrix acts from the left on spinors and from 
the right on cospinors, e.g., we resolve the Dirac operators as 

D'lp'^ = i>]^j_^^ + rml)\ Dip = -^^i)-^ + rml). 

If one strictly sticks to such convention, spinor indices can be safely suppressed, as we have 
already done in the main body of the paper and as we will often do in the remainder of this 
appendix. 

To conclude this subsection, we point out a few useful identities between the objects we have 
previously introduced. Starting from the gamma matrices, the product of an odd number of 
them has a vanishing trace. At the same time, if we consider an even number 

Tr 7^7j, = Ag^^, Tr -f^^iulalp = '^{g^iugap - gfiagup + gi^pgua) 
Tr^yall3\l[-yl5]lel^ = 4(5'[(w/pfi'/3] [7S'5]e + 5'[a£fi'/3] [557]^ + g[a5g0\-yge^) , 
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where [ ] here denotes idempotent antisymmetrisation. Furthermore, 



7^^7^ = 4/4, 7'^7«7^ = -27", 7''7V7/. = 45^^/4, 

^/.^«^/3y^^ ^ _2y/^", 7''7"7'^7^7Sm = 2(7S"7'^7^ + 7^7^7"7*). (36) 
The last equahties we shall need are 

7°C'a/3 = C'a/37^ = ^Ra/Sl'^ , [Ca/3, 77] = Ra/SfyyY , Tr Capl-ylS = -'^Ra/SfS, C^^^ = 0, 

where the equalities not involving a trace can be proved by combining the symmetry properties 
of the Riemann tensor with the anticommutation relations of the 7-matrices, e.g., 

CaPl^ = \RaPprYl^l^ = \{RapPr + Rar pph' 1^ = ^i?a/3pr (7^"7'' + 7''7^7") 

= \RaPpr {-I'l'"!^ + 2/^7'' - 2/^" - 7''7"7^ + 25''"7'') 

= \Rapl' - ICo^fiY 
^ CaPl^ = ^Raf37^. 

A. 2 On the calculus of bispinor-tensors 

The notion of bispinor-tensors heuristically boils down to consider objects which contemporary 

transform as spinor-tensors at two spacetime points. In a more sound language, they are sections 
of an outer tensor product VMMWN of two vector bundles VM, WN respectively over M and 
N. VM lEl WN is nothing but a vector bundle over M x N with, calling V and W the typical 
fibres of VM and WN, V CS) W as a typical fibre. Such a construction may seem awkward, 
but, in case M = A^, it is indeed more fundamental than the familiar tensor product bundle 
VM (8) WM, the latter being constructed out of VM M WM by pulling back via the map 
M 3 X {x,x) e M X M. 

For simplicity we will choose to collect all possible (bi)spinor-tensorial objects under the 
name of (bi)tensor, except in special case where we want to stress the character of the involved 
vector spaces. The bitensors occurring in this work are all defined only on a convex normal 
neighbourhood, since we need a unique geodesic to connect the two points our bitensors depend 
on. We use unprimed indices to indicate components stemming from tensorial properties at 
X and primed indices for those rooted in such properties at y. Furthermore, we shall use the 
bracket notation introduced by Synge to denote coincidence point limits of bitensors, namely, 

[R{x,y)] = lim B{x,y), 
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where B is some smooth bitensor, such that the hmit is weh defined. 

Let us now recall the bitensors used in this work and examine their properties. We will only 
mention the basic points while we refer to the works of DeWitt and Brehme, Fulling, Christensen, 
and to the review by Poisson for further, more exhaustive, details |DeWBr6n[lFil89llChr76[[R^ . 

As a starting point, we consider the halved squared geodesic distance cr{x, y) taken with sign, 
sometimes also called Synge 's world function. Even if the geodesic distance itself might not be 
globally smooth, it is such on geodesically convex normal neighbourhoods (provided smoothness 
of the metric) and it furthermore fulfils cr-^a''^ = 2a, an identity which can be either explicitly 
computed or derived from geometric considerations. In the following, we will, as it is customary, 
drop the semicolon when indicating covariant derivatives of a. The aforementioned equation 
together with [a] = [a^] = and [a^u] = Ofiu, two identities arising out of the defining properties 
of the geodesic distance, completely suffice to determine a, as well as all the properties we need, 
namely, the coinciding point limits of its higher derivatives. These can be obtained by means of 
an inductive procedure; as an example, in the case of [c^j^p], one differentiates fXpO"'^ = 2a three 
times and then takes the coinciding point limit. Together with the already known relations, 
one obtains [cr^j^p] = 0. At a level of fourth derivative, a new feature enters the fray, namely, 
one gets a linear combination of three coinciding fourth derivatives, though with different index 
orders. To relate those, one has to commute derivatives to rearrange the indices in the looked-for 
fashion, and this ultimately leads to the appearance of Riemann curvature tensors, i.e.. 

We stress that the discussion of these few identities is indeed much more valuable than just 
yielding the stated results since a potential reader is now able to calculate coinciding point 
limits both of derivatives of arbitrarily high order and of any bitensor; this holds true provided 
he is given the limits of lower order derivatives, an equation relating them to the higher ones, as 
well as the information of appropriate curvature tensors. We would like to remark at this point 
that, since one is ultimately interested in the coinciding point limits of certain bitensors most 
of the time, the in between computational steps often only require the knowledge of coinciding 
point limits of hierarchically lower objects, in contrast to having the necessity to know their full 
form. 

The next interesting bitensor is that of parallel transport along a geodesic, an object de- 
pending both on the underlying vector bundle and on the considered linear connection. We will 
denote the parallel transport relating the tangent spaces at x and y as g^,, while the one relating 
spinors at those same points is denoted as I^,. With them at hand, parallel transporting a 
spinor-tensor T = T^^^ ^ along the geodesic connecting y to x amounts to the following 
identity 

rpAn _ J A „tl rpB'u' 

^ — ^B'9u'^ ' 

and a similar rule applies to higher spinor-tensors. One can reverse the role of x and y, introduc- 
ing the inverses of the above two parallel transports, say gii and b • ^ practical ground, 
the construction of these two quantities boils down to finding a solution of the following partial 
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differential equations: 

ff^^'X = 4;.^" = and [gt,]=9i, [4] = I^, 

being /4 the 4x4 identity matrix. These identities, together with the properties of a and the 
inductive procedure described at the beginning of this paragraph, allow us to explicitly compute 
the derivatives of the parallel transports, the lowest ones being 

We shall henceforth suppress spinor indices, taking care to follow the afore described conventions, 
and, to conclude the section, we would like to point out the special parallel transport properties 
of both a and the gamma matrices. For the former we have, due to its geometric meaning, 

whereas, for the latter, being covariantly constant, we have 

In this paper, we need to cope with the coinciding point limits of bitensors differentiated at 
both X and y. The first, and maybe obvious, related statement is that derivatives at different 
points commute, so that we can always rearrange derivative indices in such a way that the 
unprimed ones are always on the left whereas the primed ones are always on the right. As a 
subsequent step, one notices that mixed coinciding point limits can be also calculated out of 
inductive paths. If one has the knowledge of the coinciding derivatives at the point x, however, 
one can extend it to those at y by means of Synge's rule: 

Lemma A.l. Let T be a smooth bitensor of arbitrary order; then its covariant derivatives 
possess the following property in the coinciding point limit (here suppressing all unessential 
indices): 

This has been proven by Synge for a exclusively, while, for the proof of an extension to 
arbitrary bitensors, one can refer to section 2.2 in Poisson |Po03] or to Christensen |Chr76] . 

A. 3 On the Hadamard recursion relations and related results 

As we have seen in section [3.31 in order to "construct" the two-point functions uj^{x,y) of a 
Hadamard state, we need to specify the distribution kernels H^{x,y) and the smooth bispinor 
W{x,y), which must satisfy 

P-^H^{x,y) G E{DM 0D*M), PyH^{x,y) G E{DM®D*M) 



52 



and 



D'.D'y [H^[x,y) + W{x,y)) = Py[H^{x,y) + W{x,y)) =Q. 



y 

From this it follows that D'^DyH^{x,y) is a smooth bispinor as well. Furthermore, due to 
proposition 13.21 there are even more differential operators, which, applied to H^{x,y), yield a 
smooth bispinor. Let us collect them all in the following: 

P, = -D'^D,, Py = -D'yDy, D'^D'y, D^Dy and - D'y = Dy - D',. (38) 

The aim of this section is to use these data to determine the Hadamard bidistributions H^{x, y) 
and to calculate the various coinciding point limits of their derivatives which are necessary for 
the proof of theorem 14.11 Following the path paved in the preceding sections, let us recall that 
the index structure of is 

uj^{x,y)=oj^{x,y)A^' E^{x) Es'iy), 

and that and W inherit this structure, and let us suppress spinor indices in the following. 

Although H^[x,y) and W{x,y) are bispinors, we recall from the main body of the paper 
that their form slavishly mimics that of the kernels specifying the two-point function in the 
theory of scalar fields, viz., 

HHx, y) = + V{x, y) In (39) 

oo 

V{x,y) = Y,yn{x,y)a{x,yr, (40) 

n=0 

oo 

W{x, y) = Y, ^n{x, y)a{x, y)", (41) 

n=0 

where (T±e{x, y) = (t(x, y) ± 2ie {T{x) — T{y)) + e^, with e > and T being a temporal function 
whose existence is guaranteed since the background is globally hyperbolic jBeSa051 IBeSa06| . As 
already commented in the main text, A is a reference distance employed to make the argument 
of the logarithm dimensionless, while the remaining objects, the so-called Hadamard coefficients 
U and V, are smooth bispinors. As we will see shortly, U as well as V depend only on the 
geometry of the underlying background and the mass, whereas W fully characterises the state, 
namely, the two-point functions of two Hadamard states differ only by a smooth function and 
such a difference is indeed encoded in W. 

To determine U, V, and W, we need to use the knowledge on the differential operators 
()38p which, once applied to H^, give smooth bispinors. To make the following formulas more 
readable, we choose to omit the regularising e - and thus the it index of -, the reference 
length A, and the dependence of the kernels on the spacetime points {x,y). That said, we can 
in principle take either of the second order differential operators listed in (j38p to recursively 
calculate U and V; we will employ Px, as this is familiar from the computations in the scalar 
case. 
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Applying to H, we obtain potentially singular terms proportional to cr^" for n = 1,2,3 
and to Ino" as well as smooth terms proportional to positive powers of a. We know, however, 
that the total result is smooth and one possible way to achieve this is to demand that the 
coefficients of the potentially singular terms are identically vanishing. Let us stress that, since 
we do a priori not know if U contains positive powers of a, the terms proportional to negative 
powers of a could in principle cancel each other to yield a smooth result. It is therefore a choice 
and not a necessity to require the coefficients of the inverse powers of a to vanish, and it is, 
furthermore, not guaranteed that the result of this procedure does not depend on the choice of 
the second order differential operator out of the possible ones listed in ()38p . The afore laid down 
line of argument does, however, not hold for the coefficients of PxH proportional to In a; since 
U and V are required to be smooth, they can not contain a logarithmic dependence on a and 
the terms proportional to In a have to vanish necessarily. 

The result of the previously described procedure are the so-called Hadamard recursive rela- 
tions, which, in the scalar case, have been studied by several authors (see for example |MoOO| ) . 
In the case of Dirac fields, there are results on the coinciding point limits of the Hadamard co- 
efficients up to Vi computed in |Chr78j : the form of the Hadamard singularity employed in this 
work is, however, a different one related to the non-rigorous DeWitt-Schwinger expansion, but 
formally, the relation between the different recursion relations arising in the two constructions 
is well-known. 

After having discussed the Hadamard recursion relations, we shall show how they arise 
explicitly. Let us thus examine the terms U/a and Vina individually. Starting with the latter, 
we have 

CO 

P^iVlna) = (P,y)lna + ^(K(n,a-2 + 4n) + 2a^K;^)f7"-i, 

n=0 

where we have employed the identity u'^cj^ = 2a. Remembering our previous discussion, we can 
now extract our first differential equation by requiring the coefficient of Ino" to vanish. Since 
this requirement has to hold independently of the differential operator chosen out of (jSSp . we 
have 

P^V = PyV = D'.D'yV = {D', - Dy)V = {D[ - D,)V = 0. (42) 
To obtain further differential equations, we need to look at the terms involving U, viz., 

fU\_ {P^U) , 2a^U>' + {Dxa-4)U 

which, combined with the a~^ coefficient coming from the series obtained out of differentiating 
Vina, leads us to the following two identities: 

PxU + 2Vo-^,a'' + (□..CT - 2)^0 = 0, (43) 
2U,^a^' + {n,a - A)U = 0, (44) 
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referring to the a^^ and cr~^ coefficients, respectively. 

Let us now focus on (j44p : one can infer that U is subject to a linear partial differential 
equation which, according to standard theorems, provides a unique solution once a suitable 
initial condition is given. The latter is usually chosen in such a way that 

[U] = h, 

and, hence, the Cauchy problem associated to the [/-bispinor strongly suggests us to hypothesise 
U to be of the form U = ul, with a smooth biscalar u satisfying [u] = 1. Plugging in this ansatz 
in (f44j) and recalling the properties of /, it holds that ul is the solution we are seeking if and 
only if u satisfies the partial differential equation 

2u.f,a^' + {n^a - 4)u = 0. 

Hence, it turns out that u fulfils the same transport equation as the coefficients of the 
Hadamard bidistribution encoding the singularity of the two-point function for a scalar quantum 
field and is thus given by the square root of the so-called VanVleck — Morette determinant. It 
would be tempting to think that a similar result and interpretation holds for V, but, alas, this is 
far from being the truth as one can realise by direct inspection of (j43p since spin curvature terms 
not proportional to the identity enter the arena via derivatives of /. The enlarged complexity 
of the Diracian Hadamard coefficient V, however, is compensated by the increased number of 
differential equations fulfilled by V ()42p . Of course, any of them is enough to determine V, but 
computations are still considerably easier if one employs all. 

To obtain differential equations for the Vn, one has to combine (|40p with ()42p . After a few 
formal manipulations, one gets to 

oo oo 

5^ {PM + 5^ (2Z11;^a^ + {lOa + 21 {I- 1)) Vi) a''' = 0, 

71=0 1 = 1 

and, if we require this identity to hold true at each order in a, to 

P,Vo + 2Vi.f,a'' + ia,a)Vi = 0, (45) 
P^Vn + 2{n + l)Vn+l■^,a^' + {{n + l)a^a + 2n{n + l)){Vn+i) = 0. Vn > 1 (46) 

At this point it is clear how to determine U, V and W explicitly: the starting point is ()44p . 
which, as we have explained above, gives us U once an initial condition has been assigned. 
Afterwards one can plug the result in (j43p in order to obtain Vq, though one needs to specify 
an initial condition. This is already included in (]43p . however, since, if we take the coincidence 
point limit of (jl3]) and recall the properties of o", we end up with 

m = -l[p.u]. 

Hence, we can now proceed iteratively, namely, we exploit ()45p to construct Vi once we have 
specified the initial condition taking the coincidence point limit, i.e., 

[Vi] = -\[P.Vo]. 
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Similarly, (P6|) grants us that the same procedure allows us to express Vn+i out of the preceding 
term Vn together with the initial condition 

= "2(n + l)(n + 2)[^^^"]- 

Let us remember that all these results can be obtained starting from [U] = and employing 
differential equations which only involve local curvature terms and the mass m. Thus, both U 
and V indeed only depend on these data and are independent from the state under consideration. 
This of course changes once we want to determine the final unknown quantity W . Starting from 
(j4ip and recalling the differential equation 

P^{H + W) = Q^ P,W = -P,H, 

it is clear how to get recursive differential equations for the Wn- This time two initial conditions, 
namely, [Wq] and [VFo;;^]) however, have to be specified by hand, which is of course not surprising 
since we expect some indeterminateness which has to be fixed by selecting a specific state. 

It seems that we finally have all ingredients necessary to calculate the sought coinciding 
point limits used in theorem 14.11 There is one potential feature of the Hadamard coefficients, 
however, which helps a lot simplifying calculations and should therefore be discussed before 
starting calculations, namely, their symmetry. Indeed, such a property has been proven in 
[MoOOj for the scalar case, but, unfortunately, a similar result does not exist for Dirac fields and 
even understanding the correct notion of "symmetry" in our framework is a rather challenging 
task. We shall leave the tantalising endeavour to prove the symmetry of the Diracian Hadamard 
coefficients for possible future work and circumvent, for the time being, this gap with more 
explicit calculations. The following lemma will turn out to be rather useful in general and in 
the context of coping with the lack of (proven) symmetry in particular: 

Lemma A. 2. Given a smooth bitensor B{x,y) and a smooth biscalar f{x,y) such that ^(^^'y^ 
is a smooth bitensor and 

[B] = [B,^.'] = [f] = iU'] = 0, as well as [□ J] ^ 0, 

it holds 

"Si _ [OB] 



Proof. We only sketch the proof here. Since B, f and B/f are smooth, their coinciding point 
limits do not depend neither on y nor on the path along which one approaches x. Thus, we 
can apply de I'Hospital's rule to our smooth bitensors restricted to arbitrary smooth curves ji, 
thereby expressing coinciding point limits of fractions as those of directional derivatives, e.g., 



'B{x,yy 




iiiyy' B{x,y).^> 




ii{yY' i2{yY' B{x,y).^,y, 


_ myY'i2{yr'B{x,y),^.,,] 


J{x,y) 




ii{yy f{x,y);fj,' 




ii{yyi2{yY'f{x,y)-^'y' 


[ii{yY'i2{yY' f{x,y);/i'u'] 
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where we assume that [iiiy)^' i2iyY' fi^iy);ij.'u'] is non- vanishing. This holds due to the hy- 
potheses of the lemma once we find two smooth curves 71, 72 such that 71^72^ = g^^ . 
Going to normal coordinates at x, it is always possible to find 71, 72 joining x and y such 
that 71 (x) = (1, 1, 1, 1) and 72(x) = (—1, 1, 1, 1); thus 71^^72'^ = rj^^ ■, where rj denotes the met- 
ric in normal coordinates, r/ = diag(— 1, 1, 1, 1). Since [ii{yY i2{yY fixjy);fi'iy'] is coordinate- 
independent, the statement of the lemma holds as a consequence due to Synge's rule. □ 

The last worthy of mention tool to perform the calculations whose results we will display 
shortly is the computer. It should be clear at this point that there are lot of recursion relations 
to solve to achieve the wished- for results. Thus, at the least as a means of backing up manual 
calculations, computer algebra systems are a valuable instrument. To this avail, we have chosen 
to work with Mathematica and the free package [Riccij . suitable for performing calculations with 
vector bundles. The codes we have used to implement the recursive procedures and coinciding 
point limits are available upon request from t.p.hack@gmx.de. 

We can now finally state the main proposition of the appendix: 

Proposition A.l. The Hadamard bidistribution H fulfils 

1. [P,H] = 6[Vi], [{P.H).,,] = 8[Vi.^], [(P.F);^'] = -8[Vi.^] + 6[Vi].^, 

2. [PyH] = 6[Vi], [{PyH).J = 8[Vi,,] - 2[Vi],„ [{PyH).y] = -8[Vi.J + 8[Vi],„ 

3. Tr[D',D'yH] = -Tr[P,H], Tr[{D',D'yH),^] = -Tr[{P^H)J + [V^],^, 
Tr[{D',D'yH),,,] = -Tr[{P,H),,,] - [Vi],„ 

4. Tr[{PyH - P.//)y]7'^7. = 2Tr[Vi],,. 



Proof. 

1. We shall employ (l42|), (03]), and dS]). These data entail 

P^H = 2Vi,,a<^ + Vi{U^a + 2) + 0{a), (47) 

and thus, taking the coinciding point limit and remembering those of a computed in the 
previous section, [PxH] = 6[Vi]. Similarly, one gets, deriving (|17|) once and performing 
the limit, [{PxH)-^ = 8[Vi;^]. By means of Synge's rule we finally have [{PxH)-^^i] = 

-8[yi;^]+6[yi];^.' 

2. We would of course like to compute PyH, but without any knowledge on the symmetries 
of the Diracian Hadamard coefficients, we have to verify the transport equations for Py, 
which otherwise would follow automatically from those for Px as it happens in the scalar 
case for the scalar Hadamard coefficients u and v [MoOOj. To wit, 

2U.^,>(T^'' + U{Uy(T - A) = I (2u.^,a''' + uiBya - 4)) = 0, 
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where the first equahty holds since the derivative of / vanishes along the geodesic connect- 
ing X and y and the second one holds since u{x, y) = u{y, x j^ and u is thus subject to 
transport equations for both Px and Py. Since PyH is smooth, we now know that 

^ . Yi . PyU + 2Vo.f,c7^ + Voinya-2) 

Zi = — = , 

a a 

must be smooth too. Alas, it does not factorise into a term only involving the scalar 
coefficients u and v times / and, up to now, we are unaware of a way to prove that it is 
identically vanishing. But we can try to compute whether it vanishes up to the derivative 
order we need for our purposes. To this end, it helps to split V into vl + V, where V 
is the non-trivial matrix part of V stemming from the spin curvature. This way one can 
separate from Yi a term which vanishes due to the transport equation for v and has to 
cope with the remainder only. Involved calculations yield 

[Y,] = [Y,.j = pyi] = [{nv^yj = o, 

and thus, employing lemma [Al2l [Zi] = [Zi-^] = 0. Consequently, 

PyH = 2Vi- gia^ + Vi{\Ilya + 2) + terms vanishing in the limit 

and 

{PyH)-^ = 2Vi;g'(7.^ + Vi;^(nj,(T + 2) + terms vanishing in the limit. 

One can now straightforwardly obtain [PyH] = 6[yi], [{PyH)-^] = 8[Vi;^] — 2[Vi];^, and 
\{PyH),^,] = m■,^^-ml^.■ 

3. Let us define 

Z2 = {Dx - D'y)H = {Dy - D'^)H. 

By direct inspection, 

D'xD'^H = -PxH - D',Z2. 

Again we know that Z2 is smooth and, alas, neither this quantity nor D'^Z2 turns out 
to be vanishing. Luckily enough, we can still extract some useful results at the level of 
traced coinciding point limits, at an order of derivatives high enough for our purposes. 
One computes 

Z. ^ -^(^--^> + iD.-D'y)U-ViD.-D'yy ^ 

, ^ (48) 
U(Dx - D')a Yo , ^ ' 

= 2 + - + ln(a)(Z). - D'y)V. 



^^The symmetry of u does not have to be proved in the same long way as that of v (see [MoOO] ). but it follows 
automatically by its explicit form u{x, y) = ^ detijj^^i (^1 J/)) "^Z \9{y)\~^ ■ 
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As already discussed, the last term vanishes identically and so does the first term on 
account of 

This leaves us with Z2 = Y2/a. Involved computations, employing (D^ — D'y)V = PxV = 
to exchange higher derivative terms with terms of lower derivative order in the appearing 
commutators with 7-matrices, yield 

[Y2] = [¥2;^] = [aY2] = 0, mY2yj = 6 m,^^] . 

After a few rearrangements and out of lemma \A.2\ one gets 

[Z2]=0, [Z2;,] = m,-f^]. 

Hence, [Z2-fj] is traceless due to the antisymmetry of the commutator. By means of formula 
([36]) . one can show per direct inspection that Tr[Vi]j^ji, = Tr[Vi]g^i, which entails that 
even D'^Z2 is traceless and, thus, 

Tr[D'^D'yH] = -Tr[P,H]. 

In order to compute Tr[{D'^DyH)-^] and Tr[{D'^D'yH).^i], let us consider that DyZ2 = 
D'x'Z'2 + -fx-ff — ^y^- Employing this as well as the previous results and tricks we have 
discussed in this proof, one obtains the following chain of identities 

Tr\{D',Z2\y\ = -Tn^[Z2,.^] = Tn''[Z2-y^] - Tn''[Z2-Ji^ 

= -Tr[{D'yZ2U = -Tr[iD',Z2).J + Tr[iPyH - P,H),,] 

1 (49) 
= -Tr[{PyH - P,H),,] = -Tr[Vi].,, 

= -Tr[iD'yZ2).,^']. 
We can finally use this last calculation to obtain 

Tr[{D',D'yH),,] = -TriiP.H),,] + [Vi],„ Tr[iD',D'yH).y] = -Tr[iP,H).y] - [V,],^. 

4. Inserting the previous results, we have [{PyH — PxH)y] = 2[Vi]-^. As already discussed, 
due to 7-matrix identities, tracing [Vi] with two 7-matrices amounts to a multiplication 
with the metric. Since the operations of trace and covariant derivation commute, we have 
Tr[Vi].„Y7'' = Tr[Vi].f, and thus 

Tr[iPyH - P,H).y]Y-/^ = 2[Vi];^. 

□ 

We would like to conclude this section by stating the last ingredient necessary for proving 
theorem 14.11 the coinciding point limit of Vi, viz., 

i^^l = + + TTT^ + 7^ + ^ ^4 + 



48 1152 480 720 720 / 48 



59 



A. 4 Conserved local curvature tensors 



The explicit form of the conserved local curvature tensors spanning the regularisation freedom 
of the expected stress-energy tensor is: 



J, 



K, 



7=T / Roti3R^^ d^g = -g^y{Rfj,yR^^^ + HR) — R-^j^^ + D-R^jy — 2Raj3R'^J^y 

/\q\ Og,_ip J I 



M 



\g\ Sg^,„ 



M 



— —-g^ivRafi-ySR'^^^^ + '^RajJ-yfiR'^^^u + ^RajsR^'ffu ~ ^Rafj.R% — 'iOR^^, + 2R.^^y. 

As already stated, in four spacetime dimensions, these are related as K^y = I^y — AJ^y via 
the generalised Gauss-Bonnet-Chern theorem [A195t ltHoVe74) . Furthermore, in this case, they 
all have a trace proportional to Di? and, thus, the linear combination I^y — 2>J^u is traceless. 



References 

[A195] L. J. Alty, "The generalized Gauss-Bonnet-Chern theorem," J. Math. Phys., 36 (1995), 
pp. 3094-3105. 

[Ar70] H. Araki, "On Quasifree States of CAR and Bogoliubov Automorphisms, " Publ. RIMS 
Kyoto Univ. 6 (1970/71), 385-442. 

[Ar99] H. Araki, "Mathematical theory of quantum fields," Oxford, UK: Univ. Pr. (1999) 236 

[BGP07] C. Bar, N. Ginoux and F. Pfaffle, "Wave Equations on Lorentzian Manifolds and 
Quantization" (2007) European Mathematical Society. 

[BaRa86] A. O. Barut, R. Raczka, "Theory of group representations and applications" (1986) 
World Scientific. 

[BeSa05] A. N. Bernal and M. Sanchez, "Smoothness of time functions and the metric 
splitting of globally hyperbolic spacetimes, " Commun. Math. Phys. 257 (2005) 43 
| arXiv:gr-qc/04011 12|. 

[BeSa06] A. N. Bernal and M. Sanchez, "Further results on the smoothability of Cauchy 
hypersurfaces and Cauchy time functions," Lett. Math. Phys. 77 (2006) 183 
|arXiv:gr-qc/0512095] . 



60 



[BoHi58] A. Borel, F. Hirzebruch, "Characteristic classes and homgeneous spaces. I. " Amer. J. 
Math. 80 (1958)458, 

[BrRo96v2] O. Bratteli and D. W. Robinson, "Operator algebras and quantum statistical me- 
chanics. Vol. 2: Equilibrium states. Models in quantum statistical mechanics, " Berlin, 
Germany: Springer (1996) 517 p, 

[BrOt86] M. R. Brown and A. C. Ottewill, "Photon propagators and the definition and approx- 
imation of renormalized stress tensors in curved space-time," Phys. Rev. D 34 (1986) 
1776. 

[BDF09] R. Brunetti, M. Duetsch and K. Fredenhagen, "Perturbative Algebraic Quantum Field 
Theory and the Renormalization Groups, " a rXiv;_0901.20 38 [math-ph] . 

[BFK95] R. Brunetti, K. Fredenhagen and M. Kohler, "The microlocal spectrum condition and 
Wick polynomials of free fields on curved spacetimes," Commun. Math. Phys. 180 (1996) 
633 |arXiv:gr-qc/9510056 |. 

[BrFrOO] R. Brunetti and K. Fredenhagen, "Microlocal analysis and interacting quantum field 
theories: Renormalization on physical backgrounds, " Commun. Math. Phys. 208 (2000) 
623 |arXiv:math-ph/99 03028 1 . 

[BrFr09] R. Brunetti and K. Fredenhagen, "Quantum Field Theory on Gurved Backgrounds, " 
[arXiv:090l720~63l [gr-qc] . 

[BFV03] R. Brunetti, K. Fredenhagen and R. Verch, "The generally covariant locality princi- 
ple: A new paradigm for local quantum physics, " Commun. Math. Phys. 237, 31 (2003) 
[arXiv:math-ph/0112041] . 

[ChFr08] B. Chihan and K. Fredenhagen, "The time slice axiom in perturbative quantum 
field theory on globally hyperbolic spacetimes, " Comm. Math. Phys. 287, (2009) 513 
[arXiv:0802.1642 [math-ph]] 

[Chr76] S. M. Christensen, "Vacuum Expectation Value Of The Stress Tensor In An Arbitrary 
Gurved Background: The Govariant Point Separation Method, " Phys. Rev. D 14, 2490 
(1976). 

[Chr78] S. M. Christensen, "Regularization, Renormalization, And Govariant Geodesic Point 
Separation", Phys. Rev. D 17 (1978) 946. 

[ChDu79] S. M. Christensen and M. J. Duff, "New Gravitational Index Theorems And Superthe- 
orems," Nucl. Phys. B 154 (1979) 301. 

[DaHo06] C. D'Antoni and S. Hohands, . "Nuclearity, local quasiequivalence and split property 
for Dirac quantum fields in curved spacetime, " Commun. Math. Phys. 261 (2006) 133 
|arXiv:math-ph/0106028] . 



61 



[Den82] N. Dencker, "On the Propagation of Polarization Sets for Systems of Real Principal 
Type" J. Funct. Anal. 46 (1982) 351, 

[DeWBrGO] B. S. DeWitt and R. W. Brehme, "Radiation damping in a gravitational field, " 
Annals Phys. 9 (1960) 220. 

[DFP08] C. Dappiaggi, K. Fredenhagen and N. Pinamonti, "Stable cosmological models driven 
by a free quantum scalar field," Phys. Rev. D 77 (2008) 104015 [arXivi080 1.2850 [gr-qc]]. 

[Dim82] J. Dimock, "Dirac Quantum Fields on a Manifold", Trans. Am. Math. Soc. 269 (1982) 
133, 

[DuHo72] J. J. Duistermaat and L. Hormander, Fourier Integral Operators II, Acta Mathematica 
128 (1972) 183, 

[FeVe02] C. J. Fewster and R. Verch, "A quantum weak energy inequality for Dirac fields in 
curved spacetime," Commun. Math. Phys. 225 (2002) 331 | arXiv:matirph/0105027] . 

[FlWa96] E. E. Flanagan and R. M. Wald, "Does backreaction enforce the averaged null energy 
condition in semiclassical gravity?, Phys. Rev. D 54 (1996) 6233 [ arXiv:gr-qc/9602052] . 

[Fo84] L. H. Ford, "Quantum Instability Of De Sitter Space-Time," Phys. Rev. D31, 710 (1985). 

[FoRo04] M. Forger and H. Romer, "Currents and the energy-momentum tensor in clas- 
sical field theory: A fresh look at an old problem, " Annals Phys. 309 (2004) 306 
[arXiv:hep-th/0307199j . 

[FNW81] S. A. Fulhng, F. J. Narcowich and R. M. Wald, "Singularity Structure Of The Two 
Point Function In Quantum Field Theory In Curved Space- Time. li, " Annals Phys. 136 
(1981) 243. 

[Fu89] S. A. Fulling, "Aspects of Quantum Field Theory in Curved Spacetime," London Math. 
Soc. Student Texts 17 (1989) 1. 

[Ge68] R. Geroch, "Spinor Structure of Space- Times in General Relativity. I" J. Math. Phys. 9 
(1968) 1739. 

[Ha92] R. Haag, "Local quantum physics: Fields, particles, algebras, " Berlin, Germany: Springer 
(1992) 356 p. (Texts and monographs in physics) 

[Ha03] B. C. Hall, "Lie Groups, Lie Algebras, and Representations" (2000) Springer. 

[Hol99] S. Hollands, "The Hadamard condition for Dirac fields and adiabatic states 
on Robertson- Walker spacetimes," Commun. Math. Phys. 216 (2001) 635 
|arXiv :gr-qc/990607 6]. 

[HoWaOl] S. Hollands and R. M. Wald, "Local Wick polynomials and time ordered products of 
quantum fields in curved spacetime," Commun. Math. Phys. 223, 289 (2001). 



62 



[HoWa02] S. Hollands and R. M. Wald, "Existence of local covariant time ordered prod- 
ucts of quantum fields in curved spacetime," Commun. Math. Phys. 231, 309 (2002). 
|arXiv:gr-qc/0111108j . 

[HoWa03] S. Hollands and R. M. Wald, "On the renormalization group in curved spacetime, " 
Commun. Math. Phys. 237 (2003) 123 |arXiv:gr-qc/0209029j . 

[HoWaOS] S. Hollands and R. M. Wald, "Conservation of the stress tensor in interact- 
ing quantum field theory in curved spacetimes, " Rev. Math. Phys. 17 (2005) 227 
|arXiv:gr-qc/0404074 1. 

[Ho90] L. Hormander, The Analysis of Linear Partial Differential Operators /Spinger (2000). 

[Ho94] L. Hormander, The Analysis of Linear Partial Differential Operators ///Spinger (1994). 

[Hus96] D. Husemoller, "Fibre Bundles" Springer- Verlag (1996) third edition, 

[JuSc02] W. Junker and E. Schrohe, "Adiabatic vacuum states on general spacetime manifolds: 
definition, construction, and physical properties" Ann. Henri Poincare 3 no. 6, 1113-1181 
(2002). 

[KaWa91] B. S. Kay and R. M. Wald, "Theorems on the Uniqueness and Thermal Proper- 
ties of Stationary, Nonsingular, Quasifree States on Space- Times with a Bifurcate Killing 
Horizon, " Phys. Rept. 207 (1991) 49. 

[KraOO] K. Kratzert, "Singularity structure of the two point function of the free Dirac field on 
a globally hyperbolic spacetime", Annalen Phys. 9 (2000) 475 |arXiv:math-ph/0003015|. 

[Ko95] M. Kohler: "The stress-energy tensor of a locally supersymmetric quantum field theory 
on a curved spacetime" Ph.D. Thesis (1995) Hamburg. 

[LaMi89] H. B. Lawson and M-L. Michelsohn Spin Geometry Princeton University Press (1989). 

[Lic64] A. Lichnerowicz "Champs spinoriels et propagateurs en relativite generale", Bull, de la 
S. M. F. 92 (1964) 11, 

[Li82] A. D. Linde, "The New Inflationary Universe Scenario," In *Cambridge 1982, Proceed- 
ings, The Very Early Universe*, 205-249 

[Li83] A. D. Linde, "Quantum Creation Of The Inflationary Universe," Lett. Nuovo Cim. 39, 
401 (1984). 

[LuRo90] C. Liiders and J. E. Roberts, "Local Quasiequivalence and Adiabatic Vacuum States" 
Comm. Math. Phys. 134, 29-63 (1990). 

[Mo99] V. Moretti, "Proof of the symmetry of the off-diagonal heat-kernel and Hadamard's ex- 
pansion coefficients in general C(infinity) Riemannian manifolds," Commun. Math. Phys. 
208 (1999) 283 |arXiv:gr-qc/9902034j . 



63 



[MoOO] V. Moretti, "Proof of the symmetry of the off-diagonal Hadamard/Seeley-deWitt's coef- 
ficients in C(infinity) Lorentzian manifolds by a 'local Wick rotation'," Commun. Math. 
Phys. 212 (2000) 165 [ arXiv:gr-qc/9908068j . 

[Mo03] V. Moretti, "Comments on the stress-energy tensor operator in curved spacetime, " Com- 
mun. Math. Phys. 232, 189 (2003). 

[NaOt84] A. H. Najmi and A. D. Ottewih, "Quantum States And The Hadamard Form. 2. 
Energy Minimization For Spin 1/2 Fields", Phys. Rev. D 30 (1984) 2573. 

[O107] H. Olbermann, "States of low energy on Roberts on- Walker spacetimes" Class. Quan- 
tum. Grav. 24, 5011-5030 (2007). 

[PaTo83] L. Parker and D. J. Toms, "Renormalization Group Analysis Of Grand Unified The- 
ories In Curved Space-Time," Phys. Rev. D29, 1584 (1984). 

[Po03] E. Poisson, "The motion of point particles in curved spacetime, " Living Rev. Rel. 7, 6 
(2004) |arX iv:gr-qc/0306052j . 

[Rad96a] M. J. Radzikowski, "Micro-Local Approach To The Hadamard Condition In Quantum 
Field Theory On Curved Space-Time," Commun. Math. Phys. 179 (1996) 529. 

[Rad96b] M. J. Radzikowski, "A Local to global singularity theorem for quantum field theory on 
curved space-time," Commun. Math. Phys. 180 (1996) 1. 

[Ricci] J. M. Lee, "Ricci - A Mathematica package for doing tensor calculations in differential 
geometry, " [www.math.washington.edu/~lee/Ricci] . 

[SaVeOl] H. Sahlmann and R. Verch, "Microlocal spectrum condition and Hadamard form for 
vector valued quantum fields in curved space-time, " Rev. Math. Phys. 13, 1203 (2001) 
[arXiv:math-ph/000802 9| . 

[Sa08] J. A. Sanders, "Aspects of locally covariant quantum field theory, " Ph.D. Thesis (2008) 
University of York, arXiv:0809.4828 [math-ph]. 

[Scha95] G. Scharf, "Finite quantum electrodynamics: The Causal approach, " Berlin, Germany: 
Springer (1995) 409 p. (Texts and monographs in physics) 

[St80] A. A. Starobinsky, "A new type of isotropic cosmological models without singularity, " 
Phys. Lett. B91, 99 (1980). 

[StrVeWo02] A. Strohmaier, R. Verch and M. Wollenberg, "Microlocal analysis of quantum fields 
on curved spacetimes: Analytic wavefront sets and Reeh-Schlieder theorems," J. Math. 
Phys. 43 (2002) 5514 |a rXiv:math-" ih70202003 1 . 

[Ta89] S. Tadaki, "Hadamard Regularization and Conformal Transformation" Progress of The- 
oretical Physics 81, 891, (1989) 



64 



[tHoVe74] G. 't Hooft and M. J. G. Veltman, "One loop divergencies in the theory of gravita- 
tion," Annales Poincare Phys. Theor. A 20 (1974) 69. 

[Ver94] R. Verch, "Local definiteness, primarity and quasiequivalence of quasifree Hadamard 
quantum states in curved space-time, " Gommun. Math. Phys. 160 (1994) 507. 

[Ver96] R. Verch: "Scaling Analysis and Ultraviolet Behaviour of Quantum Field Theories in 
Curved Spacetimes" Ph.D. Thesis (1996) Hamburg Universitat. 

[VerOl] R. Verch, "A spin- statistics theorem for quantum fields on curved spacetime man- 
ifolds in a generally covariant framework, " Commun. Math. Phys. 223, 261 (2001) 
[arXiv:math-ph/0102035| . 

[Vi85] A. Vilenkin, "Classical And Quantum Cosmology Of The Starobinsky Inflationary 
Model," Phys. Rev. D32, 2511 (1985). 

[Wa77] R. M. Wald, "The Back Reaction Effect In Particle Creation In Curved Space- Time, " 
Commun. Math. Phys. 54 (1977) 1. 

[Wa78] R. M. Wald, "Trace Anomaly Of A Conformally Invariant Quantum Field In Curved 
Space-Time," Phys. Rev. D 17 (1978) 1477. 

[Wa79] R. M. Wald, "On The Euclidean Approach To Quantum Field Theory In Curved Space- 
Time," Gommun. Math. Phys. 70 (1979) 221. 

[Wa84] R. M. Wald, "Ceneral Relativity" (1984) Chicago University Press. 

[Wa95] R. M. Wald, "Quantum field theory in curved space-time and black hole thermodynam- 
ics," Chicago, USA: Univ. Pr. (1994) 205 p 



65 



